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ABSTRACT 

We performed multiwavelength observations of the young planetary nebula (PN) Ml-11 and obtained its 
elemental abundances, dust mass, and the evolutionary status of the central star. The AKARIflRC, VLT/VISIR, 
and Spitzer/IRS spectra show features due to carbon-rich dust, such as the 3.3, 8.6, and 11.3 /mi features due 
to polycyclic aromatic hydrocarbons (PAHs), a smooth continuum attributable to amorphous carbon, and the 
broad 1 1.5 and 30 /xm features often ascribed to SiC and MgS, respectively. We also report the presence of an 
unidentified broad feature at 16-22 /im, similar to the feature found in Magellanic Cloud PNe with either C-rich 
or O-rich gas-phase compositions. We identify for the first time in Ml-11 spectral lines at 8.5 (blended with 
PAH), 17.3, and 18.9 /im that we attribute to the Ceo fullerene. This identification is strengthened by the fact that 
other Galactic PNe in which fullerenes are detected, have similar central stars, similar gas-phase abundances, 
and a similar dust composition to Ml-11. The weak radiation field due to the relatively cool central stars in 
these PNe may provide favorable conditions for fullerenes to survive in the circumstellar medium. Using the 
photo-ionization code CLOUDY, combined with a modified blackbody, we have fitted the ~0. 1-90 /tm spectral 
energy distribution and determined the dust mass in the nebula to be ~3.5 x 10~ 4 M Q . Our chemical abundance 
analysis and SED model suggest that Ml-11 is perhaps a C-rich PN with C/O ratio in the gas-phase of +0.19 
dex, and that it evolved from a 1-1.5 Mq star. 

Subject headings: ISM: planetary nebulae: individual (Ml-11), ISM: abundances 
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1. INTRODUCTION 

The most stable fullerene is Ceo (Draine 2011). Fullerenes 
together with other carbon dust such as graphite are expected 
to be important components of the interstellar medium (ISM), 
because they contribute to interstellar extinction. For exam- 
ple, Dopita et al. (1997) argued the possibility that the deep 
2200 A absorption feature in the low-excitation planetary neb- 
ula (PN) SMP LMC8 might be by a surface charge slop res- 
onance on Cgo- The graphite grains have a spectral peak 
around this wavelength due to it— >ir* electron excitations. 
Fullerene and polycyclic aromatic hydrocarbons (PAHs) re- 
semble graphite; therefore such grains have strong electronic 
transitions around 2200 A (Draine 201 1). The investigation of 
circumstellar carbon grains such as Ceo would be important to 
understand ISM evolution more deeply. 

Recent Spitzer/IRS studies show that fullerenes Ceo are de- 
tected in several young PNe, a proto-PN, and two post asymp- 
totic giant branch stars (post AGB) in the Milky Way (Cami 
et al. 2010; Garcfa-Hernandez et al. 2010; Zhang & Kwok 
2011; Gielen et al. 2011) and also in a handful young PNe 
in the Magellanic Clouds (Garcfa-Hemandez et al. 2011a; 
Bernard-Salas et al. 2012). The detection cases are increas- 
ing, however the excitation and formation mechanism is still 
unclear (e.g., Cami et al. 201 1; Garcfa-Hernandez 2012). The 
detections generally suggest that fullerenes can survive or be 
observed in a C-rich environment with a weak radiation field. 

Ml-11 (PN G232.8-04.7) is a good sample in C 60 forma- 
tion in a circumstellar environment, because most of the Ceo 
PNe have cool central stars (^30,000 K) and C-rich nebulae, 
and they are very young (~ 1000 yr). Indeed, Ml-11 has the 



cool central star (29 300 K, Phillips 2003), and the nebula is 
relatively compact, with a size of <6" in diameter in Ha emis- 
sion, and it appears to be a very young object (^1000 yr after 
leaving the AGB phase, this paper). Henry et al. (2010) mea- 
sured the C abundance in Ml-11 using recombination lines 
(RLs) and the O abundance with collisionally excited lines 
(CELs, number density ratio of C/O = 79.4). To date, only 
upper limits to the intensities of the CELs C III] 1906 and 
1909 A are obtained (Kingsburgh & Barlow 1994), and thus 
the gas phase C/O ratio of Ml-1 1 derived from the same type 
of emission line is still unknown. However, Ml-11 is also 
known to be a dust-rich PN, showing predominantly C-rich 
dust. Silicon carbide (SiC) and amorphous silicate features 
are seen in the IRAS/LRS data (Zhang & Kwok 1990). A sig- 
nificant near-infrared excess suggests the presence of hot dust 
(Phillips & Ramos-Larios 2005; Zhang & Kwok 1990), al- 
though there is also a contribution from the 3.3 /im emission 
feature due to PAHs in the near infrared (Allen et al. 1982). 
Longer wavelength emission features at 6.2, 6.9, 7.7, and 8.6 
/tm due to PAHs are reported by Cohen et al. (1986). A very 
weak and tentative feature is seen around 6.9-7.0 /im (see Co- 
hen et al. 1986). Indeed, there are other Ceo transitions at 
7.0 and 8.5 /im but Cohen et al. (1986) show some evidence 
(although the quality of their data is not satisfactory) for a pos- 
sible detection also in Ml-11. Therefore, Ml-11 is perhaps a 
C-rich PN. 

To confirm whether C6o in Ml-11 is real and also obtain 
insights about C6o formation, we need to investigate the phys- 
ical and chemical properties of the dust and ionized nebula 
and the nature of the central star, then we need to compare 
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the derived quantities with those in C6o PNe. To characterize 
Ml-1 1, we obtained continuous data from the UV to the far- 
infrared using several instruments, and we comprehensively 
investigated this PN. 

In this paper, we discuss C^) in Ml-11 based on the infor- 
mation of the dust and gas compositions and the evolutionary 
status of the central star. In section 2, we describe the UV- 
to mid-infrared (mid-IR) spectroscopic data from the Interna- 
tional Ultraviolet Explorer {WE), Subara/HDS, OAO/ISLE, 
AKARI/IRC, VLT/VISIR, and SpitzerfIRS, as well as narrow- 
band imaging obtained with WFPC2 on the Hubble Space 
Telescope (HST). The description in Section 2.6 includes a 
list of dust features seen in Ml-1 1, specifically, we report the 
discovery of weak features at 8.5, 17.3 and 18.9 fim lines, 
which are attributed to the Ceo fullerenes. The derivation of 
the ionic and elemental abundances in the ionized nebula is 
given in Section 3. In Section 4, we discuss the observed C 
and O gas abundances and compare them with the predictions 
from nucleosynthesis models for AGB stars. Using the photo- 
ionization code CLOUDY (Ferland et al. 1998), we fitted the 
spectral energy distribution (SED) and determined the dust 
mass and the evolutionary status of the central star. Section 4 
also includes a discussion on the formation of Cm, and a com- 
parison of the physical properties of Ml-1 1 to those of other 
Galactic PNe that exhibit fullerene features. A summary and 
future prospects are given in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Subaru/HDS observations 

We obtained optical spectra of Ml-11 using the High- 
Dispersion Spectrograph (HDS; Noguchi et al. 2002) attached 
to one of the two Nasmyth foci of the 8.2-m Subaru telescope, 
on October 6, 2008 (program ID: S08B-110, PI: M.Otsuka) 
and October 18, 2005 (PI: A.Tajitsu). The spectra were taken 
in two wavelength ranges: 3600-5400 A (the blue spectra, 
taken in 2008) and 4600-7500 A (the red spectra, taken in 
2005). 

When we obtained the blue spectra, an atmospheric disper- 
sion corrector (ADC) was used to minimize the differential at- 
mospheric dispersion over the broad wavelength region. We 
used a slit width of l."2 (0.6 mm) and a 2x2 on-chip bin- 
ning. We set the slit length to be 8" (4.0 mm), which fitted 
the nebula well and allowed us to directly subtract sky back- 
ground from the object frames. The slit position angle (PA.) 
was -225°. The CCD sampling pitch along the slit length 
projected on the sky was ~ 0."276 per binned pixel. The re- 
solving power reached around R >33 000, which is derived 
from the mean of the full width at half maximum (FWHM) 
of narrow Th-Ar and night sky lines. The total exposure 
time was 600 sec (=300 sec x2 frames). For the flux cali- 
bration, blaze function correction, and airmass correction, we 
observed G192B2B as a standard star. 

For the red spectra, we used the red image de-rotator and set 
it to PA.=90°. We set the slit width to 0."6 and the slit length 
to 7" and selected a 1 (wavelength dispersion) x 2 (spatial di- 
rection) on-chip binning. The resulting spectral resolution R 
is >65 000. We used an exposure time of 300 sec and ob- 
served G192B2B as a standard star. 

For both sets of observations, we took several bias, instru- 
mental flat lamp, and Th-Ar comparison lamp frames. We 
are interested in detecting weak C, N, and O recombination 
lines. The peak intensities of these lines are typically — 10- 
20% higher than the local continuum, and therefore a high 



signal-to-noise ratio (S/N) of the continuum is necessary. The 
resulting S/N after subtraction of the sky background was 
found to range from -5 at -3700 A to -30 at -5200 A in 
the blue spectrum, and from -5 at -4800 A to 15 at —6700 
A in the red one. 

Data reduction of the Subaru/HDS spectra and analysis 
of the emission lines was done with the long-slit reduction 
package noao . twodspec available in IRAHJ, and was per- 
formed in the same manner as described by Otsuka et al. 
(2010). When measuring the fluxes of the emission-lines, we 
assumed that the line profiles were Gaussian and we applied 
a multiple Gaussian fitting technique. 

The line fluxes were de-reddened using 
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where /(A) and F(X) are the de-reddened and the observed 
fluxes at A, respectively, and /(A) is the interstellar extinc- 
tion parameter at A, from the reddening law of Cardelli et al. 
(1989) with Ry = 3.1. The interstellar reddening correction 
was performed using the reddening coefficient c(H(3), at H/3. 
We compared the observed Balmer line ratios of H7 (blue 
spectrum) or Ha (red spectrum) with H/3 to the theoretical 
ratio computed by Storey & Hummer (1995) assuming the 
electron temperature T e = 10 4 K and the electron density « e 
= 10 4 cm -3 assuming that the nebula is optically thick in Ly- 
a (Case B of Baker & Menzel 1938). We derived c(H/3) = 
1.677±0.008 for the blue and 1.218±0.017 for the red spec- 
tra. 

The flux scaling was performed using all emission lines de- 
tected in the overlap region between the blue and the red spec- 
tra. The de-reddened relative intensities /(A) detected in both 
spectra are consistent within 10% of each other. The com- 
bined de-reddened spectrum is presented in Fig. Q] and the 
detected lines are listed in the Appendix (Table floY We have 
detected over 160 emission lines, thus exceeding the number 
of detections by Henry et al. (2010), who report more than 
70 lines in the 3700-9600 A spectra. Our measurements of 
the line intensities /(A) are in agreement with the results from 
Henry et al. (2010) within a —14 % error. 

Specifically, we detected C II, N II, and O II RLs, and highly 
excited lines due to He II, C III, and N III. These high excita- 
tion lines show a relatively broad FWHM (—0.8-1.6 A) com- 
pared to typical nebular lines (—0.2-0.5 A). It is possible that 
the He II, C III, and N III lines are not of nebular origin but of 
stellar origin, because the effective temperature of the central 
star (29 300 K; Phillips 2003) is not high enough for species 
with an ionization potential (IP) >40 eV to exist in the nebula. 
For example, we did not detect any nebular lines from species 
with an IP >40 eV, such as [Ne III] AA3876/3967 (IP>41 eV) 
and [Ar IV] A471 1/40 (> 40.7 eV). The IPs of He II, C III, and 
N III are 54.4, 47.9, and 47.5 eV, respectively. 

For our analysis, we also used the emission-line fluxes in 
the 7700-9300 A range measured by Henry et al. (2010) 
scaled in such a way that the shorter wavelength part of their 
spectrum matches our HDS observations. 

2.2. HST/WFPC2 archive data 
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omy (AURA), Inc., under a cooperative agreement with the National Science 
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FIG. 1. — The scaled and de-reddened HDS spectrum of Ml-11. The wave- 
length is shifted to the rest wavelength in air. 

We downloaded archival //5T/WFPC2 photometry in the 
F656N (6564 A/28 A) filter (PL: R.Sahai; PID: 8345), which 
traces the Ha emission. We reduced the photometric data us- 
ing the standard HST pipeline with MultiDrizzle, and present 
the drizzled Ml-11 F656N image in Fig. [2] The plate scale is 
0.025" pixel" 1 . The image shows that Ml-11 is an elongated 
nebula; the dimensions of the bright rim are —0.8" along 
PA.=-27° and -0.5" along P.A.=+63°. 

2.2.1. The total Ha and Hf3 fluxes 

The de-reddening formula (eq. (1)) requires the total H/3 
flux over the entire nebula to obtain the line fluxes in the 
AKARI and Spitzer spectra. This can be derived from the 
WFPC2 Ha image. We find that the total flux in the F656N 
filter is 1.77(-11) ± 3.48(— 13) erg s" 1 cm" 2 integrated over 
the entire PN (we will use the notation X(-Y) for X x 10" y , 
hereafter), where we assume that the uncertainty corresponds 
to the standard deviation of the background. Using the HDS 
red spectrum and the transmission curve of the F656N filter, 
we estimated —12.5 % of the total measured flux to be due to 
a local continuum and the [N II] 6548 A line flux. Thus, we 
estimated the solo Ha line flux to be 1.55(-1 1) erg s" 1 cm" 2 . 
Using the observed F(Ha)/F(H/3) ratio (6.57) in the HDS red 
spectrum, we derive the total logF(H/3) to be -1 1.629 erg s" 1 
cm" 2 , which is comparable to logF(H/3)=-l 1.84 erg s" 1 cm" 2 
measured by Cahn et al. (1992). 

2.2.2. The hydrogen density profile 

In Fig. [2a), we present the radial profile from A toward B 
indicated in Fig. [2] Based on this radial profile, we exam- 
ined the ionized hydrogen density n(H + ) as a function of the 
distance from the central star R. When we restrict the inte- 
gration to the optional portion of the nebula, the de-reddened 
observed Ha flux using c(H/3)=1.218 and the reddening law 



of Cardelli et al. (1989) with R v 
is given by: 



: 3.1, //(Ha) in erg s cm 
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where D is the distance to Ml-11 from us; j(Ha) is the emis- 
sion coefficient; and e is the filling factor which is the fraction 
of the nebular volume filled by ionized gas. V; is the volume 
of the of the nebula. In Case B, 
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FIG. 2.— HST/WFPC2 F656N image of Ml-1 1. rotated according to elon- 
gation of the nebula. In the inner 2" x2" box, the gray scale is adjusted to 
show the brightness of the bright rim and the central star. The radial profiles 
in Ha and «(H + ) from A (Dec. relative position=-0.5") towards B (-2.7") 
are presented in Fig.f3}a) and (b), respectively. 

with T e = 5400 K derived from the Balmer jump (See section 
3.2), D = 2.1 kpc (Tajitsu & Tamura 1998), and assuming n e 
~ n(H + ), n(H + ) can be written as a function of the distance 
from the central star R in cm as 



«(H + ) ~ 3.74 x 10 1 



/,(Ha) 
eR 
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The resulting n(H + ) profile from A towards B with different 
values for e is presented in Fig. 0b). If we assume that the 
ionized gas is concentrated within R=1.0" and that the den- 
sity has a constant value of 10 5 cm" 3 obtained from Balmer 
decrements, then we find that e is around 0.2. We used these 
n(H + ) profiles with different values for e in the SED modeling 
(See section 4.2). 

2.3. ISLE Observations 

We obtained / and Ks band medium-resolution (/?— 2500) 
spectra using the near-infrared imager and spectrograph ISLE 
(Yanagisawa et al. 2006, 2008) attached to the Cassegrain 
focus of the 1.88-m telescope at the Okayama Astrophysical 
Observatory. The observations were done in ISLE engineer- 
ing time in March 2008 (Ks) and January 2010 (/). The de- 
tector of ISLE is a lKx IK HgCdTe HAWAII array. We used 
a science grade detector for the /-band observations and an 
engineering grade detector for the Ks-band observations. The 
entrance slit width was 1" for both sets of observations. We 
fixed the PA. at 90° during the observations. The sampling 
pitches in wavelength were — 1.68x 10~ 4 and — 3.4x 10" 4 /j,m 
pixel" 1 in the / and Ks spectra, respectively, while the sam- 
pling pitch in the space direction was 0"25 pixel" 1 for both 
spectra. We observed standard stars HIP35132 (A0V) and 
HIP35180 (A1V) for the /-band and HIP31900 (F0V) for the 
Ks-band spectra at different airmasses to calibrate the flux lev- 
els, and correct for telluric absorption and airmass extinction. 
We observed Ml-11 in a series of 120 sec exposures in both 
observing modes. The total exposure times were 3600 sec 
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FIG. 3. — The radial profiles from A to B (Fig. [2} in Ha (upper) and «(H + ) 
(lower) with different filling factors e. 



for the /-band spectra and 3480 sec for the Ks-band spec- 
tra, respectively. Dark frames with the same science expo- 
sure time, Ar and Xe lamp frames, and on- and off-dome 
flat frames were also taken. For further wavelength calibra- 
tion and distortion correction, OH-lines recorded in the object 
frames were used. The data reduction was performed in a 
standard manner using IRAF. The interstellar reddening cor- 
rected spectra are presented in Figj4}a) and (b). For this cor- 
rection, we adopted c(H/3) applied in the HDS red spectrum. 
The resulting S/Ns are >40 in the /-band and >30 in Ks-band 
spectra at the continuum level. 

We detected more than 50 lines in these spectra, including 
a series of vibration-rotation excited lines of molecular hy- 
drogen (H2), as listed in the Appendix (Table The line 
fluxes were normalized such that /(Pa/3)=100 in the /-band 
and 7(Br7)=100 in the Ks-band spectra. 

Fig|5]shows the spatial profiles of B1-7 and H2 1-0 S(l)/2-l 
S(l) lines in the /T-band spectrum. In the both spectra, the H2 
lines are easily distinguished with other ionic lines by their 
spatial spread upto ~ 16" in the diameter. The ratio of H2 1-0 
S(l)/2-l S(l) is a traditional shock indicator (e.g, Hora & Lat- 
ter 1994, 1996; Kelly & Hrivnak 2005). The ratio in Ml-1 1 is 
~4.5 at the center of the nebula, and it decreases up to ^1.0 
outside of the ionized region. The ratio of H2 1-0 S(l)/2-l 
S(l) = 4.5 along the optical nebula (<6" in diameter) indi- 
cates a mix of UV and shock excitation (this is usual in proto- 
PNe, e.g., Kelly & Hrivnak 2005). Furthermore, the detection 
of a series of H2 lines with upper vibrational level (v > 3) in 
/-band indicates that these lines are excited by fluorescence 
through the absorption of UV photons from the central star in 
photodissociation region (PDR). 

The excitation diagram of H2 lines in the entire slit of the 
spectra is shown in Fig|6] An ortho-to-para ratio of 3 is as- 
sumed. It clearly shows that the vibrational excitation temper- 
ature (T v n,) exceeds the rotation excitation temperature (T mt ), 
indicating fluorescence emission (cf. Shupe et al. 1998 for 
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Wavelength (rest) ([im) 
FIG. 4.— ThelSLE/to-band spectra of Ml-1 1. 

BD+30° 3639; Hora & Latter 1994 for M2-9). The difference 
of the excitation temperature between different rotation levels 
(T rof ~ 2150 K for v = 1 and / ro ,-1000 K for v > 2) indicates 
H2 lines are collisionally (shock) excited in a part (the center) 
of the nebula. However, it is evident that the UV excitation 
in PDR is still dominant for the most part of H? emission in 
Ml-11. 

Any pure rotational lines in Spitzer/IRS are not detected 
(see Table|2]). 

The lines at 1.15 and 1.19 /im are identified with [P II] (in 
the 3 Pi-'Z)2 and 3 / > 2- 1 £ , 2 transitions, respectively), represent- 
ing the discovery of these lines in Ml-1 1. Adopting the tran- 
sition probabilities of Mendoza & Zeippen (1982), the colli- 
sional impacts of Tayal (2004a) and the level energy listed 
in Atomic Line List v2.05bl^| the expected [P II] 7(1.19 
^m)//(1.15 fim) ratio is 2.63 in / e =10000 K and « e =5xl0 4 
cm" 3 . The observed line ratio (2.79±0.44) agrees well with 
the theoretical value, which confirms the identification of the 
[P II] 1.15/1.19 fim lines. Our measurement of the 7([Kr III] 
2.19 fim)/I(Bij) ratio of 3.44±0.23 also agrees with Sterling 
& Dinerstein (2008; 3.22±0.26). 

2.4. IUE archival data 

The N 2+ abundance can be estimated from the N III] A1750 
line, present in archival IUE spectra which we retrieved from 
the Multi-mission Archive at the STScI (MAST). We col- 
lected low-resolution IUE spectra taken by the Short Wave- 
length Prime (S WP) and Long Wavelength Prime (LWP) cam- 
eras (file ID: SWP25846, LWP05896, and LWP05897), all of 
which were made using the large aperture (10.3 x23 arcsec 2 ). 
In these spectra, we identified the He II 1640 A and N III] 
1750 A lines. We determined that c(H/3)=0.67±0. 12, by com- 
paring the theoretical ratio of He II 7(A1640)/(A4686)=6.56 to 
the observed value, in the case of Z e =10 4 K and « e = 10 4 cm -3 
as given by Storey & Hummer (1995). The interstellar ex- 
tinction correction was made using equation (1). The flux 
measurements of the detected lines along with the normalized 

8 see http://www.pa.uky.edu/~peter/newpage/ 
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FIG. 6. — Molecular hydrogen excitation diagram from the full /- (open cir- 
cles) and K-band (filled circles) spectra. Shown are the upper state vibration- 
rotation level populations relative to that in the V = 1, / = 3 level plotted 
against the energy of the upper state in Kelvin (K). g is the statistical weight. 
An ortho-to-para ratio of 3 is assumed. The points within vibrational levels 
fall on separate lines, as expected fluorescent-excited emission. Linear fits 
to the data on each vibrational level are plotted with the derived rotational 
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values are listed in Table Q] While we did not detect the C III] 
AA1906/09 A lines in the SWP and LWP spectra, Kingsburgh 
& Barlow (1994) show 7(C m]AA1906/09) to be 6 with an 
uncertainty greater than a factor of 2. 

2.5. AKARVIRC archival data 

We analyzed the 2.5-5.5 /im prism spectra of Ml-1 1 taken 
with the Infrared Camera spectrograph (IRC; Onaka et al. 
2007) on board of the AKARI satellite (Murakami et al. 2007). 
The data were obtained as part of a mission program, PN- 
SPEC (data ID: 3460037, PI: T.Onaka), on April 11, 2009. 
The used observing window was l'x 1'. For the data reduc- 



TABLE 1 

Lines detected in the IUE spectra. 



(A) 


Ion 


/(A) 


F(\f 
(erg s 1 cm 2 ) 




/(A) 
(/(H/3)=100) 


1640 
1750 


Hell 
NlII] 


1.177 
1.154 


5.03(-14) ± 9.99(- 
1.51(-13) ± 1.52(- 


15) 
14) 


0.51 ±0.20 
1.48 ±0.51 



a logF(H/3) = -11.63 erg s~' cm 4 . 



tion, we used the IRC Spectroscopy Toolkit for the Phase 3 
data version. Fig.|7]shows the IRC spectrum with a local dust 
continuum subtracted. The S/N is >30 for the dust contin- 
uum. Several prominent lines are visible, and their central 
wavelengths are indicated by dotted lines. The line fluxes 
are listed in Tabl^D For the IRC spectra, we derived that 
c(H/3)=1.40±0.03 by comparing the observed intensity ratios 
of H 1 4-5 (Bra 4.051 pm), 4-6 (Br/3 2.625 /im), and 5-7 (Pf/3 
4.653 /^m) to H/3 and the theoretical values of Storey & Hum- 
mer (1995) for the case of 10 4 K and 10 4 cm -3 . To correct for 
interstellar reddening, we used the ratio of the extinction at 
each wavelength to the B—V color excess, A\IE(B- V), given 
by Fluks et al. (1994), in combination with the correlation 
between H/3 and the color excess, c(H/3)=l . 47 E(B -V), from 
Seaton (1979). 

In the AKARI spectra we also found an emission band at 
3.2-3.6 /xm, which may be due to aromatic and aliphatic hy- 
drocarbon species. A similar feature is seen in PN NGC7027 
and proto-PN (PPN) IRAS2 1282+5050, the latter of which 
has a [WC1 1] central star. ISO/SWS archival spectra of NGC 
7027 and IRAS 21282+5050 are shown forreference in Fig.H 
The resonance at 3.3 /im is attributed to vibrational transi- 
tions in polycyclic aromatic hydrocarbons (PAHs; e.g. Draine 
2011). We also recognize the 6.2, 7.7, and 11.3 /im reso- 
nances due to PAHs in the spectra of Ml-1 1. In particular, 
the 11.3 pm C-H out-of -plane bending mode is seen in the 
SpitzerflRS spectrum (See next section). Cohen et al. (1986) 
already detected the 6.2, 6.7, 7.7, and 8.6 /tm resonances in 
Ml-11. The 3.3 /tm band profile of Ml-11 is similar to the 
ones seen in NGC 7027 and IRAS 2182+5050, and thus we 
assume that the emission seen in Ml-11 is also due to PAHs. 
Indeed, the 3.3 /im feature was first detected by Allen et al. 
(1982), who also measured its flux. Using the theoretical in- 
tensity ratio of H I 7(5-9) to 7(5-11) = 1.86 in the case of 
r e =10 4 K and n e =10 4 cm -3 , we removed the contribution from 
H I «=5-9 to the 3.3 /im feature and estimated the flux due to 
PAHs 7(PAH 3.3 /im) to be 7.63(-12) erg s _1 cm -2 , which is 
about twice as large as the measurement by Allen et al. (1982; 
3.1(-12)ergs" 1 cm" 2 ). 

2.6. Spitzer^S and VLT/VISIR archival data 

Ml-11 was observed by Spitzer on November 10, 2006 
with the 9.9-19.6 /tm (SH) and 18.7-37.2 /im (LH) modes 
on the Infrared Spectrograph (IRS, Houck et al. 2004), as 
part of program ID 30430 (PI: H. Dinerstein; AORKEY: 
19903232). We downloaded the archival spectral images, and 
after masking bad pixels using IRSCLEAN, we extracted the 
one-dimensional spectra using SPICE. The S/N is >30 for the 
dust continuum. 

In the spectrum of Cohen et al. (1986), very weak and tenta- 
tive features are seen around 7.0 and 8.5 /tm that may arise at 
least partially from fullerene C60 with possible blending from 
the PAH 8.5 /tm band. However, the quality of their data is 
insufficient to be confident of these features. To check the 
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FIG. 7. — Comparison of the AKARIflRC spectrum of Ml-11 (upper 
line), the archival /SO/SWS spectra of the PN NGC7027 and the PPN 
IRAS21282+5050 (middle & lower). The identified lines are indicated by 
the the broken lines. The IDs are indicated by lower case letters; a: Br/3 2.63 
/im; b: H I 5-1 1 2.86 /im; c: PAH 3.29 /im C-H stretch + H I 5-9 3.29 /im; d: 
PAH 3.38/3.40 /im asymmetric CH 3 , CH 2 stretch; e: PAH 3.46 pm lone C-H 
stretch; f: PAH 3.49/3.51 fj,m symmetric CH 3 , CH 2 stretch; g: PAH 3.56^im 
aldehydes C-H stretch; h: H I 6-19,20 3.63 /im; j: Bra 4.05 /im; k: H I 4- 
5,6-14 4.05 fim; 1: He I 3-5 4.30 /im; m: H I 6-12 4.38 pm; n: Pf/3 4.65 iim 
(Reference of the wavelength of the hydrocarbon lines: Kwok 2007). 



TABLE 2 

Detected lines in the AKARIimC, Spitzer/IRS, and VLT/VISIR 

OBSERVATIONS. 





Ion 


/(A) 




F(Xf 




/(A) 


(pm) 






(erg s cm 


) 


(7(H/3)=100) 


2.63 


H I (Br/3) 


-0.955 


2.09(- 


-12) ± 2.32( 


-13) 


4.1 ±0.5 


2.86 


Hi 5-11 


-0.962 


6.12(- 


-13) ± 8.41( 


-14) 


1.2 ±0.2 


3.29 


PAH 


-0.971 


8.81(- 


-12) ± 2.55( 


-13) 


16.2 ± 1.0 




+ H I 5-9 












3.41 


PAH 


-0.973 


1.47(- 


-12) ± 2.98( 


-13) 


2.7 ± 0.6 


3.49 


PAH 


-0.974 


2.04(- 


-12) ± 3.77( 


-13) 


3.7 ± 0.7 


3.56 


PAH 


-0.975 


6.92(- 


-13) ± 3.47( 


-13) 


1.3 ±0.6 


3.63 


H I 6-19 


-0.976 


1.59(- 


-12) ± 4.08( 


-13) 


2.9 ±0.8 




+ H I 6-20 












3.73 


HI 5-8 


-0.977 


1.06(- 


-12) ± 4.58( 


-13) 


1.9 ±0.8 


3.82 


H I 6-16 


-0.978 


1.19(- 


-12) ± 8.38( 


-13) 


2.1 ± 1.5 


3.91 


H I 6-15 


-0.979 


6.25(- 


-13) ± 5.40( 


-13) 


1.1 ± 1.0 


4.05 


H I (Bra) 


-0.980 


4.94(- 


-12) ± 7.94( 


-14) 


8.9 ± 0.5 




+ H I 6-14 












4.30 


He I 3-5 


-0.982 


6.97(- 


-13) ± 6.52( 


-14) 


3.1±1.4 


4.38 


H I 6-12 


-0.982 


5.72(- 


-13) ± 7.07( 


-14) 


1.0 ±0.1 


4.65 


H I (Pf/3) 


-0.984 


8.99(- 


-13) ± 1.17( 


-13) 


1.6 ±0.2 




+ H I 6-1 1 












8.5 


C 60 * 


-0.970 


3.66(- 


-12) ± 4.59( 


-13) 


6.8 ± 0.9 


8.6 


PAH 1 " 


-0.970 


5.10(- 


-12) ± 5.07( 


-13) 


9.4 ± 1.0 


8.99 


[Arm] 


-0.959 


1.60(- 


-12) ± 1.75( 


-13) 


3.1 ±0.4 


12.37 


Hi 6-7,8-11 


-0.980 


9.68(- 


-13) ± 4.24( 


-14) 


1.0 ±0.1 


12.71 


He I 


-0.982 


4.85(- 


-13) ± 5.35( 


-14) 


0.5 ±0.1 


12.81 


[Nell] 


-0.983 


2.90(- 


-11)±4.35( 


-13) 


30.9 ± 1.5 


17.3 


C 6 o rf 


-0.981 


4.61(- 


-12) ± 4.78( 


-13) 


4.9 ± 0.6 


18.71 


[S III] 


-0.981 


1.31(- 


-12) ± 7.27( 


-14) 


1.4 ±0.1 


18.9 


C60 f 


-0.981 


6.88(- 


-12) ± 2.55( 


-13) 


7.4 ± 0.4 


33.16 


He I 


-0.993 


1.47(- 


-12) ± 2.57( 


-13) 


1.5 ± 0.3 


35.83 


He I 


-0.993 


3.89(- 


-12) ± 2.85( 


-13) 


4.0 ± 0.3 


Note. 


— For 


interstellar 


reddening correction, 


we used 



c(H/3)=1.403±0.025 for the AKARI/IRC and 1.629±0.020 for the 
Spitzer/IRC spectra, respectively. 

a logf (H/3) = -1 1 .63 erg s" 1 cm- 2 . 
b The FWHM is 0.17 /tm. 

The sum of two Gaussian components representing the PAH 8.6 fj,m. See 
text in detail. 

d The FWHM is 0.51 ±0.04 ^m. 
c The FWHM is 0.36±0.03 fim. 
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FIG. 8. — The Spitzer spectra of Ml- 1 1 (black line) andMl-12 (gray line). 
(inner box) The line-profiles of C6o 17.33 and 18.94 um. The line-profiles of 
C 60 at 8.5 /im in Ml-11 and Ml-12 are present in FigfTol 

presence of the C60 8.5 /im feature, we downloaded the 7.7- 
13.3 pm archival spectral data obtained using the VLT spec- 
trometer and imager for the mid-infrared ( VISIR) at ESO VLT 
UT3 (ID: 084.D-0868A; PI: E.Lagadec). We reduced the raw 
data using ESO gasgano. To compare the VISIR data with 
the Spitzer spectra of other C6o PNe and also to combine it 
with the Ml-ll's Spitzer spectrum, we degraded the original 
VISIR spectral resolving power of ^400 down to 90 by using 
a Gaussian convolution technique. The S/N of the convoluted 
VISIR spectrum is >70. 

The IRS and VISIR combined 7.7-37.2 pm spectrum which 
is shown in Fig. [8] reveals the solid-state/molecular features 
and atomic lines on the dust continuum thermal emission. The 
detected lines are listed in Table|2] The emission around 8.5 
/im is the complex of Ceo 8.5 pm and PAH 8.6 /im. We mea- 
sured the flux density of Ceo 8.5 /im by using multiple Gaus- 
sian fitting to separate its flux from the PAH 8.6 /im feature 
(See section 2.6.3). We performed the interstellar reddening 
correction in a similar way as described in Section 2.5. By 
comparing the observed intensity ratios H I F (n=6-l)l F(Wf3) 
& F(n =8-ll)/F(H/3) to the theoretical values of Storey & 
Hummer (1995) for the case B assumption in r e =10 4 K and 
« e =10 4 cm" 3 , we derived that c(H/3)=1.63±0.02. 

2.6.1. Broad spectral features at 10-13 and 16-22 pm 

Ml-11 appears to have C-rich dust, as evidenced by the 
presence of a broad 10-13 /im feature, which is usually at- 
tributed to SiC in the literature. This feature (centered around 
11.3 pm) is seen on top of a featureless continuum, presum- 
ably due to amorphous carbon (AC). The PAH features around 
10-11 pm are also visible in the spectrum. A second broad 
feature is seen around approximately 16-22 /im, similar in 
appearance to broad features reported in several PNe in the 
Magellanic Clouds (Stanghellini et al. 2007; Bernard-Salas et 
al. 2009; Garcfa-Hernandez et al. 201 la; Garcfa-Hernandez 
et al. 2012) 

While Stanghellini et al. (2007) associate this emission fea- 
ture with carbon-rich dust, Bernard-Salas et al. (2009) show 
that the 16-22 pm broad feature differs from the PAH plateau 
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TABLE 3 

C AND ABUNDANCES DERIVED FROM CELS, FOR MAGELLANIC 

Cloud PNe showing the 16-22 fiM feature. 



Nebula 


C 


O" 


C/O 


Ref. 


SMC1 


8.11 


8.26 


0.71 


(1) 


SMC6 


8.03 


8.22 


0.65 


(1) 


LMC8 


7.93 


8.26 


0.47 


(2) 


LMC25 


8.29 


8.17 


1.32 


(3),(4) 


LMC48 


8.40 


8.24 


1.45 


(3),(4) 


LMC85 


8.74 


8.40 


2.19 


(2) 



REFERENCES. — (1) Idiart et al. (2007); (2) Dopita et al. (1997); (3) 
Stanghellini et al. (2005); (4) Leisy & Dennefeld (2006) 
a The number density relative to the hydrogen is defined as log H=12. 



around 16-20 pm (Van Kerckhoven et al. 2000) and the 21 
/mi feature sometimes seen in carbon-rich post-AGB stars 
(e.g. Volk et al. 201 1) and PNe (Hony et al. 2001), and is more 
similar to the 18 /mi amorphous silicate feature, thus assign- 
ing a oxygen-rich carrier for this feature in objects. Bernard- 
Salas et al. (2009) imply that if the silicate identification for 
this feature is correct, the Magellanic Cloud PNe often show 
a dual dust chemistry, although earlier in the same study they 
state that not a single source in their sample shows a mixed 
chemistry, thus undermining the silicate identification of the 
16-22 /im bump. Indeed, in a comprehensive study of a large 
sample of Galactic and Magellanic Cloud PNe, Stanghellini 
et al. (2012) conclude that none of the LMC PNe considered 
show a dual chemistry, lending more credibility to the idea 
that the 16-22 /im feature is carried by a carbon-based mate- 
rial. 

Considering the gas chemistry adds to the confusion. Table 
[3] shows the C and O gas-phase abundances derived from the 
CELs for six of the PNe with 16-22 /im emission features dis- 
cussed by Bernard-Salas et al. (2009) and Garcfa-Hernandez 
et al. (201 la). Three of the six are O-rich in their gas-phase 
material, while the remaining three MC PNe actually show 
a C-rich chemistry. Since the C abundances of SMC1 and 
SMC6 have relatively large uncertainty (0.25 dex), these two 
SMC PNe could be C-rich. Even with this information, it re- 
mains unclear whether the 16-22 /imbump is due to an O-rich 
or a C-rich carrier. Thus we do not assign any identification 
to this feature, and do not include fitting the 16-22 /im bump 
in our analysis. 

2.6.2. 30 fim broad feature 

The carrier of the 30 /im feature remains somewhat of a 
mystery. While MgS has been proposed and often used as the 
carrier of this feature (Hony et al. 2003), recent work has cast 
doubt on the identification with MgS (see e.g., Zhang et al. 
2009; Garcfa-Hernandez et al. 2010; Zhang & Kwok 2011). 
Nevertheless, we consider MgS as a possible dust component 
in Ml-1 1 to explain the 30 /im feature in the spectrum. 

Other potential carriers, in particular hydrogenated amor- 
phous carbon (HAC) (see Grishko et al. 2001; Hony et al. 
2003), should also be considered. In the SED model forproto- 
PN HD 56126, Hony et al. (2003) showed that broad emission 
features around 7-9 /im and 10-13 /im are partly due to HACs, 
making the identification of the 30 /im band with HACs (Gr- 
ishko et al. 2001), a possibility. We examined whether HAC 
can be the main contributor to the 10-13 and 30 /im broad 
features by applying a modified blackbody model to the IRS 
spectrum, including both MgS and HACs. 

The observed flux density F\ due to thermal emission from 



TABLE 4 

Adopted optical constants for each model dust component. 



Dust species 


Data source 


SiC 


Pegourie (1988) 


amorphous carbon (AC) 


Rouleau & Martin (1991) 


MgS 


Begemann et al. (1994) 


HAC 


Hony et al. (2003) 



dust grains is given by 

F\ = y^ j (J^irajpiD 2 ^ m di Q\^B\(T d ,i), (5) 

where a,- is the grain radius of component i, p, is the dust den- 
sity, nid i is the dust mass, Q\j is the absorption efficiency, 
and B\(Td,i) is the Planck function for a dust temperature Tdj. 
We adopt a distance D = 2.1 kpc (Tajitsu & Tamura 1998). 
In Table H] we list the optical constants that we use for each 
of the dust species considered. For MgS, we used the optical 
constants of nearly pure MgS, e.g. Mgo.9Feo.1S, from Bege- 
mann et al. (1994). For the HACs, we adopt the Q\ of HAC 
in the case of H/(H+C)=0.3 calculated by Hony et al. (2003). 
In two models, we consider different compositions, consist- 
ing of combinations of PAHs, amorphous carbon (AC), SiC, 
MgS, and HAC (Fig.©. We considered 9.9-37.2 /im Spitzer 
spectrum and AKARI FIS 65/90 /im photometry data, except 
for the 16-22 /im broad band. We assume spherically shaped 
grains with a radius of a=0.5 /im for AC, SiC, and HACs, and 
we excluded PAHs. To reproduce the broad 30 /im emission 
using MgS, as discussed in Hony et al. (2003), we considered 
a continuous distribution of ellipsoids (CDE, e.g., Bohren & 
Huffman 1983; Fabian etal. 2001; Min et al. 2003) and calcu- 
lated Q x of CDE MgS using Eq. 18 given by Min et al. (2003). 
To simplify, we assume that the value of each ellipsoid MgS 
grain is ~47ra 3 /3, where a is 0.5 /im. 

The results of the modified blackbody fitting is shown in 
Fig. [9] and the derived Tjj, nidj, and mass fraction for each 
dust component are summarized in Table|5] Model 1 with AC, 
SiC, and MgS can explain the observed spectrum reasonably 
well, while the model 2 with HAC instead of MgS predicts an 
unseen broad emission feature around 20 /im. However, we 
should keep in mind that it is extremely difficult to charac- 
terize carbon compounds with a mixed aromatic and aliphatic 
content - such as HACs - in the laboratory because these opti- 
cal constants are strongly variable for different chemical and 
physical conditions. Jones (2012) presents HAC theoretical 
models that show the extreme variability of HAC spectra de- 
pending on parameters such as hydrogen-content, grain size, 
etc. In any case, the 20 /im feature is much weaker than the 
30 /im feature (see Ghrisko et al. 2001) and one could detect 
the 30 /im features while not detecting the 20 /im feature. At 
present, therefore, we do not completely rule out that HACs 
as a carrier of the 30 /im broad feature. 

In this paper, therefore, we assume that SiC and MgS are 
the main contributors to the 10-13 and 30 /im broad features, 
respectively. In the remainder of our analysis, we will only 
consider PAHs, SiC, AC, and MgS to model the dust emission 
in the SED. 

2.6.3. S.J, 77.5, and 18.9 \im emission due to Cm fullerenes 

In the VISIR and Spitzer combined spectrum of Ml-1 1, we 
see infrared features at 8.5 (although blended with the PAH 
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FIG. 9.— The observed SpitzerflRS and VLT/VISIR spectrum of Ml-11 
(gray lines) and the predicted SED from modified blackbody fitting (thick 
lines). Each component is indicated by a thin line. AC stands for amorphous 
carbon. 



TABLE 5 

Results from the modified blackbody fitting to the 
Spitzer/IRS SPECTRUM. 



Models 


Dust 


T d 


m d 




Comp. 


(K) 


(M e ) 


model 1 


SiC 


160 


2.85(-5) 


(FigUJi) 


AC 


120 


2.97(-A) 




MgS 


220 


3.17(-6) 


model 2 


SiC 


170 


1.75(-5) 


(Fig|9p) 


AC 


120 


1.77(-4) 




HAC 


80 


0.11 



8.6 /mi, see below), 17.3, and 18.9 /im, most likely due to the 
fullerene C6o- 

In our own Milky Way Galaxy, these C6o infrared fea- 
tures were recently detected in five PNe, including Tel, Ml- 
12, Ml -20, K3-54, and Ml -60 (Cami et al. 2010; Garcfa- 
Hernandez et al. 2010; Garcfa-Hernandez et al. 2012), a 
C-rich PPN (IPvASO 1005+79 10; Zhang & Kwok 2011), and 
two O-richpost-AGB stars (IRAS06338+5333 and HD52961; 
Gielen et al. 2011). In Figj8] we show the IRS spectrum 
of Ml-12 in comparison to Ml-11. The C 60 17.3 and 18.9 
/im features seem to be present in Ml-11 although they are 
much weaker than those in Ml-12. The spectrum of Ml-12 
shows spectral features due to PAHs, SiC, AC, MgS, and a 
very weak 16-22 /im feature, which resembles that of Ml-1 1. 
The 16-22 /im feature is not seen in Tel, Ml-20, K3-54, and 
Ml -60. Some other PNe with fullerenes in the Magellanic 
Clouds exhibit the 16-22 um feature (Garcfa-Hernandez et al. 
2011a; Garcfa-Hernandez et al. 2012). In the inset of FigH] 
the ~ 1 6-20 /im spectra of M 1 - 1 1 and Ml-12 with a local dust 
continuum subtracted, are shown. The wavelength positions 
of the intensity peak and the line widths of C6o 17.3 and 18.9 
/im lines in Ml-11 are almost coincident with those in Ml- 
12. In addition, the complex line around 8.5-8.6 /im - which 
turns out to be a blend of PAH 8.6 /im and Ceo 8.5 /im - in 
Ml-11 is very similar to that in Ml-12, as shown in Fig. [TOl 
These comparison with Ml-12 supports the identification of 
this feature as a C60+PAH blend. 

To quantify the excitation temperature and the total number 
of C60, we need to separate flux of C60 8.5 /im from the PAH 
8.6 /im band. Cami et al. (2010) reported that the FWHM of 
C60 8.5 /im in the PN Tel is 0.15 /im by Gaussian fitting. Tel 
shows strong C60 lines but very weak PAH bands. Assuming 
that the FWHM of C 60 8.5 /im is -0. 1 5 /im in M 1 - 1 1 and M 1 - 
12, we fit the broad line at 8.6 /im with multiple Gaussians. 



The line at 8.6 /im in Ml-11 can be represented by three plus 
one Gaussian components, as shown in FigfTOh. We assume 
that the profile of the PAH 8.6/im line is represented by the 
sum of two Gaussians at the peak wavelengths of ~8.7 and 
~8.8 /im. The FWHM of the PAH 8.6 /im represented by the 
sum of these two Gaussian is 0.2 /im, which is consistent with 
NGC7027 (0.23 /im in ISO/SWS spectrum shown in Fig.Q. 
The FWHM of the C60 8.5 /im line indicated by the blue line 
is 0.17 /im. The FWHM of the C 60 and PAH complex is 0.3 
/Ltm. The resultant Gaussian fitting for Ml-12 is presented in 
FigfTOb. The PAH 8.6 /im is represented by the sum of two 
Gaussians at ^8.6 and ^8.7 /im. The FWHMs of the C60 
and the PAHs for Ml-12 are as same as those of Ml-1 1. The 
measured fluxes of the Ceo 8.5 /im and PAH 8.6 /im lines in 
Ml-1 1 are listed in Table |f] 

The excitation temperature and the number of C60 were de- 
rived by creating a vibration excitation diagram as shown in 
FigHOj;. We followed the method of Cami et al. (2010). N u is 
the number of C60 molecules in the upper vibrational levels. 
N u is written by 



4irI(C 6Q )D 2 A 



(6) 



A he' 

where 7(C6o) is the fluxes of C60 lines in erg s" 1 cm" 2 , D is 
the distance to Ml-11 (2.1 kpc; Tajitsu & Tamura 1998), A 
are the transition probabilities (4.2, 1.1, 1.9 s -1 for C60 8.5, 
17.3, and 18.9 /im, respectively, from Garcfa-Hernandez et 
al. 201 lb), h is Planck's constant, and c is the speed of light. 
The vibrational degeneracy is given by g u . In thermal equilib- 
rium, the Boltzmann equation relates the N„ to the excitation 
temperature T ex! : 



N u oc g u exp(-Eu/ kT ext ). 



(7) 



Eu and k are the energy of the excited level and the Boltz- 
mann constant, respectively. We confirmed that our measured 
excitation temperature, T ex , of 338±9 K in Tc 1 using Eqs. 
© and © is consistent with Cami et al. (2010; 332 K). Ac- 
cordingly, we obtained T eK , of 399±36 K and the NiC^o) of 
4.57±1.23(+46). The total mass of C 60 m C6 o is 2.75(-8)M Q . 

Our estimated T eKI , N(C(,o), and mca) in Ml-12 are 345±35 
K, 5.30±1.23(+46), and 3.18(-8)M Q , respectively, adopting 
D of 3.9 kpc (Tajitsu & Tamura 1998). Garcfa-Hernandez et 
al. (2010) measured the T ext in Ml-12 of 546 K. Their mea- 
sured FWHM of Ceo 8.5 /im is 0.237 /im. The T ext discrep- 
ancy between theirs is due to the differences in the measured 
line flux of this line. 

3. RESULTS 
3.1. CELs plasma diagnostic 

We determined the electron temperatures T e and densities 
n e using 1 1 diagnostic CELs, and listed the results in Table [7] 

For [Nil] A5755 and [On]AA7320/30, we subtracted the 
recombination contamination from both lines using 



7 R ([NII]A5755) 



: 3.19 



T 



0.33 



N 
H+ 



2+ 



and 



7 R ([OII]AA7320/30) 

W) 



= 9.36 



T 

104 



0.44 



o 

H+ 



2+ 



(8) 



(9) 



given by Liu et al. (2000). Adopting N 2+ and 2+ io nic a bun- 
dances derived from N II and O II lines (see Section [3~4i i. we 



The detection of C 60 in the PN Ml-1 1 



9 



0.5 




8.2 8.4 8.6 8.8 9 9.2 8.2 8.4 8.6 8.8 9 9.2 
Wavelength (rest) (jum) Wavelength (rest) (jum) 




102 



100 



500 



1000 



1500 
E u /k (K) 



2000 



FIG. 10. — {panels (a) and (b» The line-profiles of C 60 at 8.5 /xm in Ml-11 
and Ml-12 (blue lines). The local dust continua are subtracted. The gray 
lines are the observations. We fit the broad line at 8.6 fim by three or four 
Gaussian. The deconvolved profiles are indicated by the black and blue lines. 
The red lines are the sum of these components. The emission-line around 9 
is [Ar III] 8.99 fim. (panel (c)) The excitation diagram for Cgo in Ml-1 1 
and Ml-12. The filled and open circles are the observed data in Ml-11 and 
Ml-12, respectively. The lines indicated represent the best fit to the data. See 
text in detail. 



— 1.4xl0 4 
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FIG. 1 1 . — n t -T c diagram. Each curve is labeled with an ID number given 
in Table [JJ The solid lines indicate diagnostic lines of the T s . The broken 
lines indicate diagnostic lines of the n e . 



TABLE 6 

Adopting n e and T e for the CEL ionic abundance 

CALCULATIONS. 



Zone 



Ions 



n e (cm 3 ) 



7e(K) 



N°,0° 

N + , + , P+, S + , Cl + 
2+ , Ne+, S 2+ , Cl 2+ , 
Ar 2+ , Fe 2+ , Kr 2 * 



1380 
36490 
51 120 



9240 
8410 
9540 



derived that I R ( [Nil] A5755) = 0.12 and I R ([0 n] AA7320/30) 
= 2.11. The recombination contamination is —2 % in [Nil] 
A5755 and —8 % in [Oil], respectively. 

The resulting n e -T e diagnostic diagram is shown in 
Fig. [TT] The solid lines indicate diagnostic lines for 
the electron temperatures, while the broken lines are 
electron density diagnostics. Since the gas in Ml-11 
has much higher density than the critical densities of 
the density sensitive lines [Oil] AA3726/29 and [S II] 
AA6717/31, the ratios of [On] /(AA3726/29)//(AA7320/30) 
and [S II] /(AA6717/31)//(AA4069/76) are density-sensitive 
rather than temperature-sensitive. The critical densities 
at r e =10000 K for [Oil] AA3726/29 are -4500 and 
-980 cm -3 , respectively, and those of [S II] AA6717/31 
are —1400 and —3600 cm -3 , respectively. The densities 
derived from the [On] AA3726/29 and [Sn] AA6717/31 
ratios might be the value for thin shell region, while 
those from the [O n] /(AA3726/29)//(AA7320/30) and [S n] 
7(AA6717/31)//(AA4069/76) would be the value for the bright 
rim. To determine n e , we adopted r e =10000 K for all density- 
diagnostic lines. r e ([Nll]) was calculated using n e =36490 
cm -3 , which is the average between « e ([Oll]) derived 
from/(AA3726/29)//(AA7320/30)(n e ([O ll]„ /fl ) hereafter) and 
« e ([Slll]). For r e ([Slll]), r e ([Olll]), and r e ([Arlll]), we 
adopted n e ([Slll]). For J e ([Ol]), we adopted n e ([Nl]). Our 
estimates for T e ([N II]), r e ([0 III]), and r e ([S III]) are in excel- 
lent agreement with those by Henry et al. (2010), who esti- 
mated 10720, 9996, and 9100 K, respectively. The estimates 
for n e and T e are summarized in Table [7] 

To calculate the CEL ionic abundances, we adopt a three 
zone model based on n & -T e diagram. The T e and n e combi- 
nations for each ion are listed in Table [6] For the N° and O 
abundances, we adopted r e ([Ol]) and n e ([N I]) (zone 1). The 
averaged n e is from n e ([0 ll]„/ fl ) and n e ([S III]); and T e ([N II]) 



is for ions with 0<IP<13.6 eV (zone 2). n e ([Slll]) and 
the averaged temperature among r e ([Olll]), r e ([Slll]), and 
T e ([Ar III]) is for ions with I.P.= 14.4-35.5 eV (zone 3). 

3.2. RL plasma diagnostics 

We calculate the He, C, N, and O abundances using RLs 
of these elements by adopting the T e derived from the Balmer 
discontinuity together with He I line ratios and the n e from 
the Balmer decrements, listed in Table|7] The Balmer discon- 
tinuity temperature T e (BJ) was determined using the method 
described by Liu et al. (2001), which we used to obtain the 
C 2+ , N 2+ , and 2+ abundances from recombination lines. 

The He I electron temperatures T e (He I) were derived from 
the ratios of He 1 7(A7281)//(A6678), 7(A7281)//(A5876), and 
7(A6678)//(A5876) assuming a constant electron density of 
10 6 cm" 3 , estimated from the Balmer decrements (see below). 
We adopted the emissivities of He I by Benjamin et al. (1999). 
The r e (He I) derived from three different line ratio combina- 
tions is 3980-6980 K. We adopted 7/ e (He I) derived from He I 
7(A7281)//(A6678) for the He + abundance calculations. The 
reason why this T e (He I) is the most reliable was discussed in 
Otsukaetal. (2010). 

The intensity ratios of the high-order Balmer lines Hn (n: 
the principal quantum number of the upper level) to a lower 
Balmer line, ej* H/3, are also sensitive to the electron den- 
sity. In Fig. [12J we plot the ratios of higher-order Balmer 
lines to H/3 compared to the theoretical values by Storey & 
Hummer (1995) for r e (BJ) and n e of 10 5 , 5xl0 5 , and 10 6 
cm -3 . The electron density in the recombination line emitting 
regions seems to be >10 5 cm. Care when dealing with this 
value is necessary, because it apparently has large scatter. 

3.3. CEL ionic abundances 
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TABLE 7 
Plasma diagnostic results. 





ID 


Diagnostic 


Value 


Result 




(1) 


[N II] (A6548/83)/(A5755) 


42.45±1.14 <i 


8410±90 


(K) 


(2) 


[O III] (A4959/5()07)/(A4363) 


182.94±20.32 


9740±330 




(3) 


[Arm](A7135)/(A5192) 


132.44±26.92 


10 060±780 




(4) 


[0 1] (A6300/63)/(A5577) 


74.62±14.31 


9240±640 




(5) 


[S III] (A9069)/(A6313) 


9.96±2.63 


8830±980 






Hei(A7281)/(A6678) 


0.20±0.01 


6890±330 






He 1 (A7281)/(A5876) 


0.052±0.002 


5920±230 






He I (A6678)/(A5876) 


0.26±0.01 


3980±160 






Average 




5600 




BJ 


hi aimer Jump 




5400±1300 


lie 


(6) 


[Nl] (A5198)/(A5200) 


1.63±0.13 


1380±350 


(cm" 3 ) 


(7) 


[0 II] (A3726)/(A3729) 


2.75±0.05 


5750±380 




(8) 


[On] (A3626/29)/(A7320/30) 


3.11±0.10* 


21 860±780 




(9) 


[S ll](A6716/31)/(A4069/76) 


0.63±0.02 


35 350±1400 




(10) 


[S ill] (A18.71^tm)/(A9069) 


0.44±0.11 


51 120±15 550 




(11) 


[Clm](A5517)/(A5537) 


0.40±0.10 


22110-63 730 






B aimer decrements 




10 3 -10° 



:l Corrected recombination contribution for [N II] A5755. 

b Corrected recombination contribution for [O II] AA7320/30. 
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FIG. 12. — Plot of the intensity ratio of the higher order Balmer lines to 
H/3 (Case B assumption) with the theoretical intensity ratios in 7/ c =5400 K 
and different n c . 

The derived ionic abundances are listed in Table [8] In the 
last line of the transition series for each ion, we present the 
adopted ionic abundance in bold face. The adopted values 
represent the line intensity weighted mean in case two or 
more lines are detected. As references, the results by Ster- 
ling & Dinerstein (2008, for Kr) and Henry et al. (2010, for 
the others) are also listed in the last column. This is the first 
time the Ne + , P + , and Fe 2+ abundances are derived for Ml- 
ll. In total, fifteen ionic abundances are determined by solv- 
ing for a >5 level atomic model, with the exception of Ne + , 
for which the abundance was calculated using a two level en- 
ergy model. We adopted the same collisional strengths and 
transition probabilities used in Otsuka et al. (2010, 201 1) ex- 
cept for Cl + , for which we adopted the transition probabilities 
from the CHIANTI atomic databasfl the collisional impacts 
of Tayal (2004b), and the level energy listed in Atomic Line 
List v2.05bl2. We subtracted the recombination contamina- 
tion in the [Nll]A5755 and [Oll]AA7320/30 lines to derive 
the N + and + abundances. 

In general, our derived abundances are comparable to the 
values of Henry et al. (2010) and Sterling & Dinerstein 

9 http://www.ukssdc.ac.uk/solar/chianti/ 



(2008). The discrepancies between our N + , + , S + abun- 
dances and the results by Henry et al. (2010) are mainly due 
to adopted electron temperature; Henry et al. adopted 10200 
K. Note that they employed three different r e s and a con- 
stant « e =20000 cm" 3 in their models. With their r e s, we 
would find N + , + , and S + abundances comparable to their 
results. Lastly, a minor discrepancy for the N + and + abun- 
dances must be due to the recombination contamination in the 
[N n]A5755 and [O n] AA7320/30 lines, respectively. 

3.4. RL ionic abundances 

Our values for the ionic abundances derived from RLs are 
listed in Table [9] In general, the Case B assumption applies 
to lines from levels having the same spin as the ground state, 
and the Case A assumption applies to lines of other multi- 
plicities. In the last of the line series of each ion, we present 
the adopted ionic abundance and the error estimated from the 
line intensity-weighted mean. Effective recombination coef- 
ficients for the lines' parent multiplet are the same as those 
used by Otsuka et al. (2010). The RL ionic abundances are 
insensitive to the electron density under <10 8 cm" 3 (Zhang 
& Liu 2003). For the ionic abundance calculations, we there- 
fore adopted the effective recombination coefficients in case 
of n e =10 4 cm" 3 for C 2+ and N 2+ . For He + and 2+ calcula- 
tions, we adopted n e =10 6 cm" 3 and 10 4 cm" 3 , respectively. 
Since He II, C III, C IV, N III, and O III appeared to be of stel- 
lar origin, we did not estimate the abundances of these ions. 

The He + abundances are determined using electron density 
insensitive six He I lines to reduce intensity enhancement by 
collisional excitation from the He 2s 3 S level. For the C 2+ 
abundances, the V6 and V17.06 lines, which have higher an- 
gular momentum as upper levels, seem to be unaffected by 
both resonance fluorescence by starlight and recombination 
from excited 2 S and 2 D terms. Comparison of the C 2+ abun- 
dances derived from V6 and V17.06 lines indicated that the 
observed C II lines would have less population enhancement 
mechanisms. Our estimated He + and C 2+ abundances are con- 
sistent with Henry et al. (2010; 3.56(-2), 4.48M)). 

We estimated the 2+ abundances using the O II lines show- 
ing the least contamination from other ionic transitions. We 
excluded the O II 4676.23 A line when determining 2+ abun- 
dance, because this line is much stronger than the other VI 
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TABLE 8 

Ionic abundances derived from CELs. 





(A or fim) 


/(A) 
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2.20(0) ± 4.08(-2) 


5.79(-6' 


ib 1 .34-1 
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ib 2.98 1 
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rvf 


1.15 
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1 1C\/fW _L A AC/ T\ 
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ib 2. 19 1 
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A I/}/ 1 1 _L c. cc/ 91 
0.j4(— JJ it J.Jj(- L) 


A HA( 71 
4. / 4^— / ; 
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It 4. JO' 


-a) 






6716.44 


4.42(-l)± 1.61(-2) 


3.63(-7) 


ib 1.72( 


-8) 
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^770 s i 

O / ju. o 1 


1 ARfAl 4- 7 ftSC 71 


A 1 S*' 71 


4- 1 


—5; 
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4.34(-7) 


ib 2.37' 
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a 7 n in 
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1 77/ 1 1 _i_ i o\ 


Q /1^/ 71 


4- 7 77/ 
it Z.ZJ' 


71 


/.D/(-/) 
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4- 7 77/ 
It Z. / J' 


71 
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Q SQf 1 71 
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ib 9.08 
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Cl + 


8578.69 
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ib 7.18' 
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ib 3.92' 
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7 QS/ 71 4- ^ 77/ 71 
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5537.60 


7.22(-2) ± 1.15(-2) 


1.55(-8) 

l.j^y— 0, 


± 3.85( 

4- \ Qd, 


-9) 

01 
-fj 




Ar 2+ 


5191.82 


1.86(-2) ± 3.76(-3) 


2.91 (-7) 


± 1.07( 


-7) 






7135.80 


2.46(0) ±4.1 8(-2) 


2.20(-7) 


±2.51( 


-7) 


3.90(-7) 




8.99 


3.06(0) ib 3.76(-l) 


4.38(-7) 


±5.85( 


-8) 










2.41(-7) 


± 3.95( 


-8) 


1.86(-7) 


Fe 2+ 


4701.53 


6.89(-2) ± 8.46(-3) 


9.31(-8) 


± 2.28( 


-8) 






4754.69 


3.86(-2) ± 8.65(-3) 


1.26(~7) 


± 3.90( 


-8) 






4769.43 


5.32(-2)±6.12(-3) 


2.08(-7) 


± 5.03( 


-8) 






4881.00 


4.93(-2)ib5.38(-3) 


3.74(-8) 


± 9.27( 


-9) 






4934.08 


1.63(-2)±3.54(-3) 


1.49(^7) 


± 4.42( 


-8) 






5270.40 


7.06(^2) ± 3.70(-3) 


8.20(-8) 
1.09(-7) 


± 1.63( 
± 2.72< 


-8) 
-8) 




Kr 2+ 


2.19 


8.84(-2) ± 7.28(^3) 3.52(-9) ± 3.87(-10) 


2.99(-9) 



NOTE. — (in fifth column) All values except for Ki^ + are from 
Henry et al, (2010). The Kr 2+ abundance is from Sterling & Diner- 
stein (2008). 

a Based on [O Il]AA3726/29 complex. 



lines. 

The abundance discrepancy factor (ADF), which is the ratio 
of RL to CEL abundances, is 19.5±8.4for0 2+ and 15.3±13.6 
for N 2+ . The large uncertainty of ADF(N 2+ ) is due to the large 
uncertainty in the CEL N 2+ abundance. ADF(0 2+ ) in Ml-11 
is comparable to that of high density PN M2-24 (namely, 17; 
Zhang & Liu 2003). Zhang & Liu (2003) also found that the 
density structure for Ml-1 1 and M2-24 is similar, with M2-24 
showing a large density contrast of « e =10 3-6 cm -3 . 



TABLE 9 

Ionic abundances derived from RLs. 



X+ 


Multi. 


Alab 
(A) 


7(A) 
[7(H/3)=100] 


X „ 1+/H + 


He + 


Vll 


5876.66 


6.16(0) ± 8.32(-2) 


4.17(-2) ± 2.93(-3) 




V14 


4471.47 


1.72(0) ±3.02(-2) 


3.37(-2) ± 1.98(-3) 




V45 


7281.35 


3.19(-1) ± 1.27(-2) 


3.60(-2) ± 2.74(-3) 




V46 


6678.16 


1.60(0) ± 2.78(-2) 


3.83(-2) ± 1.97(-3) 




V48 


4921.93 


4.64(-l) ± 7.25(-3) 


3.38(-2) ± 1.99(-3) 




V51 


4387.93 


2.38(-l) ± 8.44(-3) 


3.78(-2) ± 2.99(-3) 
3.93(-2) ± 2.58(-3) 


c- + 


V2 


6578.01 


4.99(-l) ± 1.25(-2) 


6.76(-^t) ± 3.11 C^) 




V3 


7231.33 


2.16(-1) ± 1.25(-2) 


5.13(-4) ± 2.26(^1) 




V3 


7236.42 


4.36(-l) ± 1.08(-2) 


5.77(^1) ± 2.52(^1) 




V3 


7237.17 


7.14(-2) ± 4.99(-3) 


8.44(^4) ± 3.74(^1) 




V6 


4267.15 


6.14(-l)±2.68(-2) 


5.31(-4) ± 2.00(^1) 




V17.06 


5342.43 


3.20(-2) ± 4.30(-3) 


C n/ A \ 1 O /II t A\ 

5.37(^1) ± 2.42(^1) 
5.59M) ± 2.31(-4) 




V3 


5666.63 


1.01(-1) ± 4.80(-2) 


5.09(-4) ± 3.61(-4) 




V3 


5677.66 


5.18(-2) ± 1.00(-2) 


5.51(-4) ± 3.09(^1) 




V5 


4630.54 


5.49(-2) ± 7.05(-3) 


2.74(^1) ± 1.46(^1) 




V19 


5017.22 


1.38(-2) ± 3.58(-3) 


5.12(-4) ± 3.07(^1) 
4.61M) ± 2.92(-4) 


o 2+ 


VI 


4641.81 


2.61(-2) ± 1.01(-3) 


1.52(-4)±4.68(-5) 




VI 


4676.23 


8.59(-2) ± 6.78(-3) 


9.22(-4) ± 1.19C-4) 




VI 


4649.14 


2.20(-2) ± 6.00(-3) 


1.17(-4) ±4.80(-5) 




VI 


4650.84 


3.10(-2) ± 7.00(-3) 


2.00(^1) ± 7.59(-5) 
1.61(-4) ± 5.85(-5) 



3.5. Elemental abundances 

The elemental abundances are estimated using an ioniza- 
tion correction factor, ICF(X), based on the IP. ICFs(X) for 
each element are listed in Table [TOl The He abundance is the 
sum of the He + and He 0+ abundances, and we allowed for the 
unseen He abundance. The C abundance is the sum of the 
C + and C 2+ abundances, and we added the unseen C + using 
ICF(C). Henry et al. (2010) used ICF(C)=0/0 2+ . Since the 
IPs of C + 2+ (11.3 and 24.4 eV) are close to those of N +2+ 
(14.5 and 29.6 eV), we instead adopted ICF(C)=N/N 2+ . The 
N abundance is the sum of N + and N 2+ . For the RL N abun- 
dance, we corrected for the unseen N + assuming (N/N 2+ )cel 
= (N/N 2+ ) RL . The O abundance is the sum of the + and 2+ 
abundances. For the RL O abundance, we assume (0 2+ /0)rl s 
= (0 2+ /0)cels- The Ne abundance is equal to the Ne + abun- 
dance. The P abundance is the sum of the P + and P 2+ abun- 
dances. We assumed that P/P + = S/S + and considered the P 2+ 
abundance. The S abundance is the sum of the S + , S 2+ , and 
S 3+ abundances. We corrected for the unseen S 3+ abundance 
using the CEL O and + abundances. We assume that the CI 
abundance is the sum of Cl + and Cl 2+ . The Ar abundance is 
the sum of the Ar + and Ar 2+ abundances, and we corrected 
for the unseen Ar + . The Fe abundance is the sum of the Fe 2+ 
and Fe 3+ abundances, and we corrected for the unseen Fe 3+ . 
The Kr abundance is the sum of the Kr + , Kr 2+ and Kr 3+ abun- 
dances, and we corrected for the unseen Kr + and Kr 3+ . 

The resulting elemental abundances are listed in Table [TTJ 
The types of emission lines used for the abundance estima- 
tions are specified in the second column. The number den- 
sities of each element relative to hydrogen are listed in the 
third column, and the subsequent two columns are the num- 
ber densities in the form of log 10 (X/H), where H is 12, and the 
relative number densities to the solar value. The last columns 
are the measurements by Sterling & Dinerstein (2008) for Kr 
and Henry et al. (2010) for the others. Our estimated abun- 
dances are in agreement with the values given by these authors 
except for C, N, and O. The C discrepancy between Henry 
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TABLE 10 

Adopted Ionization Correction Factors (ICFs). 



X 


Type 


ICF(X) 


X/H 


He 


RL 


s/s- + 


ICF(He)He + 


C 


RL 


(N/N 2+ ) CEL 


ICF(C)C 2+ 


N 


CEL 


1 


ICF(N)(N + +N 2+ ) 




RL 


(N/N 2+ ) CEL 


ICF(N)N 2+ 


O 


CEL 


1 


ICF(0)(0 + +0 2+ ) 




RL 


(0/0 2+ ) C EL 


ICF(0)0 2+ 


Ne 


CEL 


1 


ICF(Ne)Ne+ 


P 


CEL 


(S/S+) 


ICF(P)P+ 


S 


CEL 


[l-(l-(0 + /0)) 3 ]- 1/3 


ICF(S)(S + +S 2+ ) 


CI 


CEL 


1 


ICF(C1)(C1+ + C1 2+ ) 


Ar 


CEL 


(N/N 2+ ) C EL 


ICF(Ar)Ar 2+ 


Fe 


CEL 


(0/0+) 


ICF(Fe)Fe 2+ 


Kr 


CEL 


C1/C1 2+ 


ICF(Kr)Kr 2+ 



TABLE 11 
Elemental abundances of M 1 - 1 1 . 



X 


Type 


X/H 


log(X/H)+12 


[X/Hf 


Others 


He 


RL 


5.65(-2)±7.36(-3) 


1075±0.06 


-0.1 5 ±0.06 


10.55 


C 


RL 


2.19(-3)±1.60(-3) 


9.34±0.32 


+0.95±0.41 


9.78 


N 


CEL 


1.22H0±1.85(-5) 


8.08±0.07 


+0.25±0.13 


7.67 




RL 


1.81(-3)±1.58(-3) 


9.26±0.38 


+1.43±0.60 




O 


CEL 


4.13M)±2.40(-5) 


8.62±0.03 


-0.07±0.06 


7.88 




RL 


8.05(-3)±3.45(-3) 


9.91±0.19 


+1.22±0.21 




Ne 


CEL 


4.76(-5)±2.98(-6) 


7.68±0.03 


-0.19±0.10 




P 


CEL 


7.51(-8)±8.59(-9) 


4.88±0.06 


-0.58±0.06 




S 


CEL 


1.43(-6)±9.34(-8) 


6.15±0.03 


-1.04±0.05 


5.96 


CI 


CEL 


2.83(-8)±5.56(-9) 


4.45±0.09 


-0.81±0.11 


4.11 


Ar 


CEL 


9.85(-7)±6.15(-7) 


5.99±0.27 


-0.81±0.11 


5.37 


Fe 


CEL 


1.13(-7)±2.97(-8) 


5.05±0.11 


-2.42±0.12 




Kr 


CEL 


6.48(-9)±2.21(-9) 


3.81±0.15 


+0.53±0.17 


4.25 



a The solar abundances are from Lodders (2003). 



et al. (2010) and ours is due to the different ICF(C). When 
we assume ICF(C)=0/0 2+ , we find that the C abundance is 
10.45±0.22. The N and O discrepancies are due to the N + 
and + discrepancies caused by the differently adopted r e 's. 

4. DISCUSSION 

4.1. C and O abundances and the C/O abundance ratio 

We attempted to estimate the carbon abundance from CELs 
by correlating ADF(C 2+ ) with ADF(0 2+ ). We found a tight 
correlation, ADF(C 2+ )= 0.997 xADF(0 2+ ), among 56 Galac- 
tic PNe from Wang & Liu (2007), Wesson et al. (2005), Liu 
et al. (2004), and Tsamis et al. (2004). When we consider 
only the 18 PNe with ADF(0 2+ )>5, there seems to be no 
correlation between ADF(C 2+ ) and ADF(0 2+ ). For these 18 
PNe ADF(C 2+ ) is 7.65±1.97. When we adopt this value for 
ADF(C 2+ ), the C 2+ abundance derived from CELs is 7.31(- 
5)±3.56(-5) and the /(C ra]A1906/09) is 33.8±16.5. The 
CEL C abundance and the CEL C/O ratio are 8.46±0.34 dex 
and -0.14±0.34 dex, respectively. The CEL C/O ratio de- 
rived from CELS agrees well with the value derived from 
RLs, which is -0.57±0.37 dex, although there is a large un- 
certainty. 

We should keep in mind that the >9 dex C abundances de- 
rived from both the CELs and the RLs are confirmed in many 
PNe and can also be theoretically explained by AGB nucle- 
osynthesis if we assume a small minimum H-envelope mass 
for the third dredge-up (~0.5 M Q ; Straniero et al. 1997 for 
solar metallicity). However, for stars with a main sequence 



TABLE 12 

Comparison of the observed abundances with the prediction 

BY AN INITIALLY 1 .5 Mq STAR WITH Z=0.004 (SEE TEXT IN DETAIL). 



He 


C 


N 


O 


Ne 


P 


Sources 


10.75 


9.34/8.46 


8.08 


8.62 


7.68 


4.88 


This work 


10.97 


8.46 


7.65 


8.23 


7.42 


4.86 


Karakas (2010) 



mass of ~l-3 Mq and an initial metallicity of Z=0.004, O 
abundances of >9 dex are not expected to occur. For exam- 
ple, we present the predicted abundances for initially 1.5 Mq 
stars with in Table fT2l In Ml-11, we adopt the observed N 
and O abundances derived from CELs, while for the C abun- 
dances we used CEL predictions and observed values for the 
RLs. Karakas (2010) adopted scaled-solar abundances as the 
initial composition. The accuracy of the predicted abundances 
is within 0.3 dex. The model predicted C and O abundances 
are close to the predicted CEL C and the observed CEL O 
abundances. The model could also explain the observed abun- 
dance of Ne and P, which are He-rich intershell products that 
are brought up to the stellar surface by the third dredge-up in 
late AGB phase. 

Comparison of the model of Karakas (2010) and our ob- 
servation suggests that the CEL O abundance would be more 
reliable in Ml-11 relative to the RL O. This would be also 
applied to the N abundance. Therefore, we consider the CEL 
CNO abundances to represent the gas-phase abundances of 
these elements in Ml-11. If the CEL C/O ratio is correct, 
Ml-11 might be a C-rich or O-rich PN. In the next section, we 
verify the C and O abundances through gas+dust SED model- 
ing. 

4.2. SED modeling 

We investigate the fractional ionization, the evolutionary 
stage of the central star, and the physical conditions of the 
ionized gas and dust grains, by constructing an SED model 
based on the photo-ionization (P-I) code CLOUDY c08.00 and 
modified blackbody fitting. Through SED modeling, we de- 
rive the ionized gas mass m g , dust mass m^, and dust temper- 
ature Td in the nebula, as well as the total luminosity L* and 
the effective stellar temperature T e ff . We verify CEL C and O 
abundances through our model. 

4.2. 1 . Model approach 

Using Cloudy and modified blackbody fitting, we at- 
tempted to fit the observed SED from ^0. 1 to 90 /im, assum- 
ing that the dust in Ml-1 1 is composed of PAHs, amorphous 
carbon (AC), SiC, and MgS grains. We excluded the broad 
16-22 /im feature in the fitting procedure. 

Since there are no optical constants of MgS available for 
the UV, we could not include MgS in the CLOUDY model- 
ing. Instead, we fitted the MgS feature using single temper- 
ature thermal emission applied to the infrared opacity. For 
MgS, we assumed CDE, while for all other dust species we 
used spherical grains. Unfortunately, CLOUDY does not al- 
low CDE grains, providing a second reason to model the MgS 
contribution separately. Indeed, at present, there are no pub- 
licly available P-I codes which can handle CDE grain shapes. 

We assume that the observed SEDs can be expressed as the 
sum of the calculated SED from the CLOUDY model and the 
modified blackbody fit to the MgS emission. First, we fitted 
the SED in the range from ^0.1 to 22.5 /im and 65-90 /im 
using Cloudy. Then, we fitted a modified blackbody to the 
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difference between the observations and the fitted SED from 
Cloudy, to account for the contribution from MgS. In both 
the Cloudy model and the MgS modified blackbody fitting, 
we adopted a distance D to Ml-1 1 of 2. 1 kpc. 

In the CLOUDY calculations, we used Tlusty's non-LTE 
theoretical atmosphere modeH with [Z/H]=-0.81, which is 
consistent with the observed [Ar/H]. McCarthy et al. (1997) 
found that the effective temperature r e ff=29 000 K, the surface 
gravity log g=3.0, and the core mass M*=0.74 M© for the cen- 
tral star, using their non-LTE model based on Keck/HIRES 
spectra. However, they adopted a distance D of 4.6 kpc. 
Guided by their results, we used a series of theoretical atmo- 
sphere models with T e ff in the range from 27 500 to 35 000 
K and log g values of 3.0, 3.25, 3.50, 3.75, and 4.0 to de- 
scribe the SED of the central star. We found that T e ff=3 1 950 
K, log g=3.5, and a total luminosity L*=4510 L© can accom- 
modate the observations well. 

For the gas-phase elemental abundances X/H, we adopted 
the observed values listed in Table QT| as a first guess and 
varied these to match observations, except for C. For C, 
we varied the C abundance to match the predicted value: 
/(C m]A1906/09)=33.8±16.5. We did not fit the He n, C in 
or N III lines, because these lines would be of stellar wind ori- 
gin. For the N abundances, we adopted the CEL value. Treat- 
ing the collisional impacts of [P II] as a function of electron 
temperature and solving the 5-level energy model for this ion 
is not included in CLOUDY c08.00, leading to an overestimate 
of the P abundance (log 1() P/H + 12 >6). We adapted the code 
to perform this analysis correctly, and obtain a better estimate 
of the P abundance. We also revised the C III], [N II], [O II, 
III], [Ne II], [S III], and [Ar III] line-calculation programs in 
Cloudy. The abundance of elements that were not observed 
was fixed at [X/H]=-0.81. Reliable dielectric recombination 
(DR) rate measurements do not exist for low stages of ion- 
ization of S at photo-ionization temperatures (CLOUDY c08 
manual), thus we adopted the scaled DR of oxygen for sulfur 
line calculations, to match the observed [S II]. 

We adopted the hydrogen density profile as presented in 
Fig. Ob). Based on the HST image shown in Fig. [2] we fixed 
the outer radius R out at 2.6" (0.026 pc) and varied the inner 
radius and filling factor e to match the observed SED. The 
best fit R in and e are -0.8" (0.008 pc) and 0.25-0.3, respec- 
tively, assuming that the gas and dust co-exist in the same 
sized nebula. 

We assumed that amorphous carbon and SiC are present in 
the form of spherical dust grains with a radius a=0.5 /an. To 
calculate the opacities for these species, we used the optical 
constants listed in Table |4] Zhang & Kwok (1990) observed a 
nearlR-excess, also shown in the SED of Fig. [13] arising from 
both the PAH and the small dust grain emission. To explain 
the excess, we also included small carbonaceous grains with 
a grain size of a=0.005 /im, as well as PAHs. We adopted the 
optical data of Desert et al. (1990), Schutte et al. (1993), and 
Bregman et al. (1989). To explain the MgS feature around 30 
/xm, we performed a modified blackbody fitting with a single 
dust temperature Tj and CDE shaped grains with a character- 
istic size of a=0.5 /im as we mentioned in Section 2.6.2. 

To verify the degree of modeling accuracy, we evaluated the 
gas emission line strengths and the broad band fluxes in the 
features of interest (Tablefl4l. including three AKARI/FIS and 
two Wide-field Infrared Survey Explorer (WISE) photometry 
bands; and we fitted the overall SED shape with the adopted 

10 http://nova.astro.umd.edu/Tlusty2002/tlusty-frames-cloudy.html 



elemental abundances. We will discuss the resulting elemen- 
tal abundance and the ionization correction factors. 

4.2.2. Results of SED modeling 

The derived physical quantities are listed in Tables Qj] and 
n~4l which can be compared with the flux levels in the 13 spec- 
tral features of interest. In Table H4] we compare the predicted 
emission line and band fluxes with the observed values, where 
/(H/3) is 100 and the intrinsic log /(H/3) is -10.22 erg s _1 
cm" 2 . Columns 4 and 5 list the observed values and the val- 
ues predicted by CLOUDY modeling. The type of emission is 
indicated in column 3. 

The predicted CEL C and O abundances are close to the 
model results by Karakas (2010). The predicted CEL C/O 
ratio (+0. 19 dex) suggests that Ml-1 1 could be a normal low- 
mass C-rich PN slowly evolving toward higher effective tem- 
peratures. The elemental abundances predicted by CLOUDY 
are comparable with the observation and are in excellent 
agreement with the model prediction by Karakas (2010). 

In Fig. [T3j we present the calculated SED from CLOUDY 
(dots), the MgS modified blackbody fitting (broken line), and 
the sum of those two components (thick line). The gray 
lines and circles represent the observations including the IUE 
spectra, 2MASS JHKs (Ramos-Larios & Phillips 2005), the 
AKARI 9/18/65/90 ^m, and the WISE 3 and 4 fim photometry. 
In the inset, we focus on the —10-40 fim part of the SED. The 
predicted SED matches the data in UV to mid-IR range well, 
except for the PAH emission, in particular the 13.6/im feature 
due to the C-H out of plane bending mode (quartet). 

The temperatures of the 0.5-fim sized amorphous carbon, 
SiC, and MgS are 91, 102 (by Cloudy), and 160-200 K 
(by blackbody fitting), respectively. The temperature of the 
0.005-/im sized amorphous carbon grains and PAHs are 177 
K and 267 K, respectively. The gas C mass within 2.6" pre- 
dicted by Cloudy is 4.68(-5) M©. Most of the C in the 
nebula is existing as the grains. The mass of pure Ceo in the 
nebula is only —0.008 % of the total C dust. 

Assuming a distance of 2.1 kpc, the initial mass of the pro- 
genitor is estimated to be 1-1.5 M© based on the location of 
Ml-11 with respect to the post-AGB H-burning evolutionary 
tracks with Z = 0.004 (Vassiliadis & Wood 1994), as shown in 
Fig. Q3] The age is > 1000 yrs after leaving the AGB phase. 

4.3. Fullerene formation in PNe 

Since the discovery of fullerenes C6o and C70 in a C-rich PN 
Tel (Cami et al. 2010), the number of detections of these lines 
are increasing (Garcfa-Hernandez et al. 2010, 2011a; Zhang 
& Kwok 201 1); however, the formation process of fullerenes 
in evolved stars and in the ISM is under debate. At present 
three explanations for fullerene formation are proposed. 

4.3.1. Destruction of HAC 

Garcfa-Hernandez et al. (2010) suggested that PAHs and 
fullerenes may be formed by the photochemical processing of 
HAC in H-rich circumstellar envelopes. In laboratory exper- 
iments, Scott et al. (1997) showed that PAHs and C 50 ,C 6 o, 
and C70 may be produced by the decomposition of HACs. 
This fullerene formation scenario seems to be supported by 
the Spitzer observations of H-poor R Coronae Borealis (RCB) 
stars (Garcfa-Hernandez et al. 2011b) In particular, the mid- 
IR spectrum of the R Coronae Borealis star V854 Cen - an 
H-rich RCB star - evolved from HACs (ISO/SWS spectrum in 
1996) to PAHs and C 60 (Spitzer/IRS in 2007). The ISO/SWS 



14 



Otsuka et al. 




Wavelength (//m) 

FIG. 13. — The fitted SED from the CLOUDY modeling (dots) and the 
modified blackbody fitting (long dash) and the resultant SED (thick black 
line). The gray circles and lines are data from WE, Subaru/HDS, OAO/ISLE 
spectra, 2MASS (Ramos-Larios & Phillips 2005), AKARIHRC spectra and 4 
bands (9/18/65/90 ^m), and SpitzerflRS spectrum. In the CLOUDY model, 
we considered PAHs, amorphous carbon (AC), and SiC. The modified black- 
body fitting is performed for MgS only. The close-up feature of the observed 
and fitted SEDs around 10-40 fim are presented in the inner box. See main 
text for a detailed description. 

1 1 1 1 1 1 1 1 1 1 1 1 1 , 1 
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1.5 M (M,=0.64 M e ) 
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2.5 - \ 1 



2 I I i i i .i i i i i I i i i i I 

5.5 5 4.5 4 

Log (T eH ) (K) 

FIG. 14. — The location of Ml- 1 1 on the 1.0 and 1.5 Mq hydrogen burning 
evolutionary tracks after AGB phase. The location of Ml- 11 is indicated by 
the filled circle. 

spectrum of this object is similar to a laboratory spectrum 
of HAC at 773 K (Lambert et al. 2001). Later, Garcfa- 
Hernandez et al. (2011a) show that seven Magellanic Cloud 
PNe have broad features around 6-9 and 10-14 /im. The cen- 
tral stars of these PNe are H-rich. The 6-9 /im broad feature 
is different from PAH features in the same wavelength range 
and and is rather similar to that of HACs. These MC PNe have 
cool central stars (3 1 300-43 300 K, meaning young PNe), and 
their UV spectra, except for SMP SMC 1 3 and LMC2, show P- 
Cygni profiles, suggesting the presence of a stellar wind. This 
suggests that the shocks may be triggering the HAC decompo- 
sition and forming both fullerenes and PAHs. More recently, 
in the strong C 60 PNe Tel, SMC 16, and LMC56, which are 
showing HAC features, Bernard-Salas et al. (2012) argued 
that the UV radiation from the central star is photochemically 
processing the HACs to produce fullerenes, by comparing be- 



TABLE 13 

The derived properties of the central star, ionized nebula, 
and dust by the sed model. 



Parameters 


Central Star 


U 


4710 L Q 


Tell 


31 830 K 


l°gS 


3.5 cm s~ 2 


[Z] 


-0.81 


M» 


~0.6 Mq 


Progenitor 


1-1.5 M 


Age 


>1000yr 


Distance 


2.1 kpc 


Nebula 


Abundances 


He:ll.ll,C:8.49,N:7.89,O:8.30, 


(logX/H+12) 


Ne:7.67,P:5.21,S:6.15,Cl:4.20, 




Ar:5.95,Fe:5.01,Others:[X/H]=-0.81 


Geometry 


Spherical 


Shell size 


R in =0.8" (0.008 pc)/R oul =2.6"(0.027 pc) 


input density profile 


See FiglJb) 


e 


0.25 


logF(H/3) 


-10.22 erg s _1 cirr 2 (de-redden) 


m g (CLOUDY) 


0.023 M Q 


Dust in Nebula 


Composition 


PAHs,AC,SiC,MgS 


Grain size 


0.5 & 0.005 fim (See text in detail) 


Grain shape 


CDE for MgS, spherical for the others 


^(Cloudy)" 


90-267 K (See text) 


r rf (BB FITTING) 1 ' 


160-200 K 


«Jj(CLOUDY) fl 


3.44H1) M e 


m d (BB FITTING)* 


4.81(-6)-9.27(-6)M 


m d (TOT.f 


3.49Hl)-3.53(-4) Mq 


m d /m g (CLOUDY) 


1.52(-2) 



NOTE. — The accuracy of log(X/H) + 12 is 0.3 dex. AC stands 
for amorphous carbon. 

a The temperature and total mass of PAHs, AC, and SiC grains 
derived by CLOUDY. 

b The temperature and mass of MgS grains derived by the modified 
blackbody fitting. 

The total mass of PAHs, AC, SiC, and MgS grains. 

tween the observed intensity ratios of C^o 7.0, 8.5, 17.3 to 
18.9 /im and the theoretical predictions (see their Figs. 5 and 
6). In Tel, the intensity peak of C 60 8.5 /im is 6400-9700 AU 
from the central star, while that of dust components is nearby 
the central star. It would be interesting that the intensity of the 
11.2 /mi feature peaks ^3000 AU away from the central star 
on the other side of of Ceo 8.5 /im (see their Fig. 3), because 
HACs seem to have a spectral peak around 11.2 /im (Jones 
2012). 

4.3.2. Destruction of PAH clusters 

Cami et al. (2010) argued that fullerenes are not formed 
from HACs, because if fullerenes formed by HAC decompo- 
sition process, fullerenes should be observed more frequently 
in objects with PAHs. Cami et al. (2011) suggested that 
fullerenes are products of the destruction of large PAHs (>60 
C atoms) by the stellar wind shock (^100 km s" 1 ). PAHs 
first lose their peripheral H-atoms in the shock and in the post 
shock gas, and then the resulting carbon clusters can assemble 
into fullerenes. 

4.3.3. Photo-chemical process of PAH clusters 

In the reflection nebula NGC 7023, Sellgren et al. (2010) 
show that the Ceo 18.9 /im emission peaks on the central star, 
while the PAH 16.4 /im emission is brightest between the 
regions of strong 18.9 /im and the H2 0-0S(l) 17 /im emis- 
sion. Berne & Tielens (2012) found that the abundance of Ceo 
and PAHs in this nebula increases and decreases, respectively, 
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TABLE 14 

The predicted relative fluxes by the Cloudy models. 



Ions 


"lab 


Tvnp 


1 V'vobs 


'v-'vcloudy 


He+ 


4387.93 A 


RL 


2.38(-l) 


2.69(-l) 




4471.47 A 


RL 


1.72(0) 


2.22(0) 




4921.93 A 


RL 


4.64(-l) 


5.84(-l) 




5876.66 A 


RL 


6.16(0) 


6.76(0) 




6678.16 A 


RL 


1.60(0) 


1.74(0) 




7281.35 A 


RL 


3.19(-1) 


4.70(-l) 


C 2+ 


1906/09 A 


CEL 


<5.03(+l) 


4.52(4-1) 


N + 


5754.64 A 


CEL 


5.87(0) 


6.92(0) 




6548.04 A 


CEL 


5.89(+l) 


6.77(41) 




6583.46 A 


CEL 


1.90(+2) 


2.00(42) 


N 2+ 


1750 A 


CEL 


1.41(0) 


1.04(0) 


o+ 


3726.03 A 


CEL 


5.54(+l) 


6.97(41) 




3728.81 A 


CEL 


1.97(+1) 


2.94(41) 




7320/30 A 


CEL 


2.58(+l) 


4.39(41) 


o 2+ 


4361.21 A 


CEL 


1.37(-1) 


1.80(-1) 




4958 91 A 


CEL 


6.44(0) 


6.07(0) 




JUUU.04 rt 


CEL 


1 85(4-1) 


1 83(41) 


Ne + 


12.81 fim 


CEL 


3.09(4-1) 


3.18(41) 


P+ 


1.15 fim 


CEL 


1.22(-1) 


1.31(-1) 




1.19 fim 


CEL 


3.41(-1) 


3.93(-l) 


S + 


4068.60 A 


CEL 


1.79(0) 


1.27(0) 




4076.35 A 


CEL 


6.34(-l) 


4.07(-l) 




6716.44 A 


CEL 


4.42(-l) 


4.56(-l) 




6730.81 A 


CEL 


1.08(0) 


9.08(~1) 


s 2+ 


6312.10 A 


CEL 


3.23(-l) 


5.50(-l) 




9068.60 A 


CEL 


3.22(0) 


3.89(0) 




18.71 /im 


CEL 


1.40(0) 


3.65(0) 


C1+ 


8578.69 A 


CEL 


1.19(-1) 


5.94(-2) 




9123.60 A 


CEL 


7.76(-2) 


1.57(-2) 


Cl 2+ 


5517.66 A 


CEL 


2.88(-2) 


3.24(-2) 




5537.60 A 


CEL 


7.22(-2) 


7.62(-2) 


Ar 2+ 


5191.82 A 


CEL 


1.86(-2) 


2.62(-2) 




7135.80 A 


CEL 


2.46(0) 


3.04(0) 




8.99 fim 


CEL 


3.06(0) 


1.91(0) 


rc 


4701 SI A 


CEL 








4754.69 A 


CEL 


3.86(-2) 


2.70(-2) 




4769.43 A 


CEL 


5.32(-2) 


2.22(-2) 




4881.00 A 


CEL 


4.93(-2) 


7.07(-2) 




5270.40 A 


CEL 


7.06(-2) 


8.85(-2) 


Bands 


^center 


AA 


/(A)obs 


7(A) C l 0Ut Jy 


n 
D 


U.4U fim 


U.aUU fim 


4.Uj(+z) 


A ] 1\ 

4.jo(4z; 


.1 


1.24 fim 


0.162 fim 


1.14(4-2) 


1.24(42) 


H 


1.66 fim 


0.251 fim 


9.21(+1) 


1.12(42) 


Ks 


2.16 fim 


0.262 ^tm 


8.36(41) 


8.01(41) 


WISE1 


3.353 fim 


0.663 fim 


1.15(4-2) 


1.32(41) 


WISE2 


4.603 


1.042 fim 


1.43(4-2) 


1.27(42) 


AKARIl 


9.22 fim 


4.104 fim 


1.62(4-3) 


1.01(43) 


I RSI 


10.30 fim 


1.000 ^tm 


2.21(4-2) 


2.18(42) 


IRS2 


11.30 fim 


1.000 Aim 


6.54(4-2) 


6.81(42) 


IRS3 


12.50 


1.000 fim 


5.14(4-2) 


5.14(42) 


IRS4 


13.50 fim 


0.600 A*m 


2.08(4-2) 


2.56(42) 


AKARI2 


66.70 fim 


20.2 fim 


4.49(4-2) 


5.04(42) 


AKARI3 


89.20 fim 


39.9 Atm 


3.70(4-2) 


3.47(42) 



when approaching the ionizing star HD200775 (r e ff— 19000 
K, e.g. Alecian et al. 2008). They proposed that fullerenes 
are products of the photo-chemical processing of large PAHs. 
Such large PAH clusters could be from HACs as described 
above. 



4.3.4. Comparison between Cm PNe 

Which explanations are more suitable to describe fullerene 
formation in Galactic PNe? To answer this question, we 
checked the properties of the central stars and the dust com- 
ponents, as summarized in Table [131 Garcfa-Hernandez et al. 



(2010) detected fullerenes in K3-54; however, the properties 
of the central star and gas-phase elemental abundances are 
unknown. Therefore, we do not list this PN. The progenitor 
masses M prop . and the evolutionary age after AGB phase in the 
fifth and sixth columns are estimated from the location on the 
H-burning tracks of Vassiliadis & Wood (1994) with Z=0.02 
for Tel and Z=0.004 for the others. Prior to our study, The C 
and Ne abundances in Ml -12 were not well known. However, 
using the measured T e {[S III]) and n e ([S II]), the intensities of 
the C II 6462 A, and the de-reddened H/3 flux given by Henry 
et al. (2010), and our measured F([Ne II] 12.8 fim) = 6.83(- 
12) erg s" 1 cm" 2 , we estimated C and Ne abundances. To 
do so, we assumed that C/H = ( Ar/Ar 2+ ) x C 2+ /H + and Ne/H 
= Ne + /H + . We derived C 2+ /H + and Ne + /H + as 2.06(-3) and 
5.96(-5), respectively. In the gas-phase abundances of all ob- 
jects listed in Table[L5] the N, O, Ne, S, CI, Ar, and Kr abun- 
dances are derived from CELs, and the He and C abundances 
from RLs. The CEL C abundances in Ml -12 and Ml -20 are 
unknown due to a lack of UV spectra. 

We found the following common properties in these PNe: 
(1) relatively cool central stars except for Ml -20, (2) young 
age, (3) broad 6-9, 10-14, and 30 fim features, (4) a strong 
stellar wind. In the UV spectra taken by IUE and FUSE, 
we found that Tel shows the P-Cygni profile of Lya 1215 A. 
There are no data with enough signal-to-noise ratio to check 
for the presence of P-Cygni Lya profiles in the other PNe. 
However, since these PNe have the Wolf-Rayet-like central 
stars ([WC10,11] and weak emission line central star, wel), 
which show broad Cm, IV lines (>60 km s" 1 at FWHM), 
there would be a strong wind for each. 

From properties (1) and (2), the short time duration of 
weak radiation from the central star would be an essential 
condition for fullerenes to survive or to be observed for a 
long time. The observational fact that all fullerene PNe 
have cool central stars and low-excitation nebulae, except for 
Ml -20, matches the detection of the C6o in a proto PN and 
two post AGB stars by Zhang & Kwok (2011) and Gielen 
et al. (2011), respectively; these sources have cool central 
stars (21500 K in IRAS0 1005+79 10; 6250 and 6000 K in 
IRAS06338+5333 and HD52961, respectively). Property (3) 
implies that fullerenes might not be formed from HACs in 
Galactic PNe. HACs reproduce 6-9, 10-14, -20, and -30 
/im broad features. If the 30 fim feature is from HACs, we 
should see the —20 fim feature, too. However, not all Galac- 
tic fullerene PNe show the —20 /im feature. As we mentioned 
in Section 2.6.2, however, the strength of 20 and 30 /im fea- 
tures largely depends on H-content, grain size, etc. At present, 
therefore, we cannot rule out that fullerenes arise from the de- 
struction of HACs. From property (4), the strong stellar wind 
would promote fullerene synthesis for young PNe and post- 
AGB stars with cool central stars. In Ml-1 1, as we mentioned 
in Section 2.3, the H2 lines within the optical diameter are 
excited by both UV fluorescence and shocks. Lumsden et al. 
(2001) observed the HK-bznd spectra of Ml -12 and Ml -20. 
However H2 1-0S(1) and 2-lS(l) lines are not detected. 

If C60 is from large PAH clusters, the spatial distribution of 
PAHs and C 6 o would be similar to that in NGC7023. If C 6 o 
is from HACs, the spatial distribution of C60 would be sim- 
ilar to that in Tel, as shown in Fig. 3 of Bernard-Salas et al. 
(2012). At present, destruction of HACs seems to be the most 
plausible scenario for fullerene formation in Galactic PNe, in- 
cluding Ml-1 1. If so, the spatial distribution of Ceo in Ml-1 1 
would be similar to that in Tel . To advance our understanding 
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TABLE 15 

Physical parameters of C 6 o detected Galactic PNe. 



Nebula 


T - 


T 


l°g g 


M a 
i" prop. 


Age 


Type ot 


Gas-phase abundances^ 


Dust 


T_ 

JC60 


Kcl. 




(K) 


(*o) 


(cm s ) 


(M@) 


(yrs) 


Central star 


(He/C/N/O/Ne/ ) 
(P/S/Cl/Ar/Kr) 




(K) 




Ml-ll 


31950 


4510 


3.5 


1-1.5 


>1000 


[WC10-11] 


10.75/9.34/8.08/8.62/7.68/ 
4.88/6.15/4.45/5.66/3.81 


PAHs,SiC?,AC,MgS? 


399 


(1),(2),(3) 


Ml-12 


33 000 


8700 


3.4 


-1.5 


~1000 


[WC10-11] 


10.54/9.40/7.55/8.11/7.78/ 
■ ■ ■ /6.25/5.01/5.40/3.83 


PAHs,SiC?,AC,MgS? 


345 


(1),(2),(3),(4),(5) 


Tel 


34700 


6900 


3.4 


~2 


<1000 


Of(H) 


> 10.78/8.56/7.56/8.42/7.80/ 
5.30/6.45/4.97/6.71/ ■■■ 


(weak)PAHs?,AC,MgS? 


339 


(1),(4),(6),(7),(8),(9) 


Ml -20 


53 000 


8300 


5.0 


1-1.5 


<3900 


wel 


10.98/8.74/7.94/8.58/8.21/ 
■■■ /6.68/4.93/5.85/ ■■■ 


PAHs?,SiC?,AC,MgS? 


425 


(1),(4),(7),(10),(11),(12) 



REFERENCES. — (1) This work; (2) Weidmann & Gamen (201 la); (3) Sterling & Dinerstein 2008; (4) Zhang & Kwok (1993); (5) Henry 
et al. (2010); (6) Cami et al. (2010); (7) Garcfa-Hernandez et al. (2010); (8) Pottasch et al. (2011); (9) Weidmann & (2011b); (10) Wang & 
Liu (2007); (11) Gorny et al. (2009); (12) Kaler & Jacoby (1991) 

a The values are estimated from the location on the H-burning tracks with Z=0.02 for Tel and Z=0.004 for the others by Vassiliadis & Wood 
(1994). 

b The number density relative to the hydrogen is defined as log H=12. The C abundance is derived from RLs except for Tel. The C in Tel 
is derived from CEL. The N, O, and Ne abundances are derived from CELs. The predicted abundances in Ml-ll by CLOUDY are listed in 
Table[[3] 



of what kinds of conditions in which fullerenes can be formed 
and survive, the high-resolution spatial observations are nec- 
essary to investigate the spatial variation of HACs, PAHs, Ceo, 
and dust components. 

5. SUMMARY AND FUTURE WORK 

We performed multiwavelength observations for the young 
PN Ml-1 1 to investigate elemental abundances and dust mass, 
and we discussed its evolutionary status from its progenitor. 
We found a large discrepancy between the RL and the CEL 
O & N abundances. The RL C/O abundance ratio is <1. If 
the RL C/O ratio represents the chemistry of the PN, the neb- 
ula is O-rich. The RL O abundance is much larger than that 
predicted by AGB nucleosynthesis models, while the CEL O 
abundance is close to that predicted by the models. The CEL 
C and O would be reliable in Ml-1 1. Therefore, we estimated 
the CEL C abundance using the relation between ADF(C 2+ ) 
and ADF(0 2+ ) among 56 Galactic PNe, and we obtained a 
CEL C/O, implying the possibility that Ml-ll could be C- 
rich. We estimated the CEL C and O abundances through the 
SED model. 

In the AKARI/IRC, VLT/VISIR, and SpitzerflRS spectra, 
3.3, 8.6, and 11.3 /m PAH bands and broad 10-13 and 30 /im 
features are visible. SiC and MgS could be the main contrib- 
utors to the broad 10-13 and 30 /im features, respectively. We 
detect broad emission from ^16-22 pm; however, it is still 
unclear whether this feature is from C-rich or O-rich dust. We 
also detect three C6o lines in the VISIR and Spitzer data of 
Ml-ll. The presence of C6o might be explained by destruc- 
tion of HACs in this PN. As shown in Table [15] fullerene- 
containing Galactic PNe seem to have central star properties, 
gas-phase abundances, and dust composition in common. Ex- 
cept for Ml -20, most of these PNe that have fullerenes seem 
to have cool central stars. The radiation field around a rela- 
tively low temperature star is intrinsically weak. Moreover, 
the surrounding shell structure around such a low tempera- 
ture star might have diluted the stellar UV flux. Such cir- 
cumstances must have played a favorable role for fullerenes 



to survive in the circumstellar medium. 

Through ^0.1-90 /mi SED modeling, we estimate the dust 
mass and the properties of the central star and the nebula. Our 
current analysis indicates that the progenitor would be a 1-1.5 
Mq star. The observed abundances and the predicted CEL C 
and C/O can be explained by an AGB nucleosynthesis model 
for 1.5 Mq progenitors with Z=0.004. Ml-ll is possibly a 
C-rich PN. 

For a future study on Ml-ll, high spatial resolution spec- 
troscopy or narrow-band imaging would be necessary to re- 
solve the spatial distribution of PAHs, Cgo, HAC, and dust 
components to check how these co-exist in the nebulae. It is 
essential to investigate the properties of the central stars, the 
excitation degree of the nebulae, and the dust composition in 
young and low-excitation PNe such as Ml-ll to understand 
the condition in which fullerenes form and survive. 
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APPENDIX 
OBSERVED LINE LISTS 

The detected lines in the Subaru/HDS and OAO/ISLE are listed in Tables 16 and 17, respectively. 



The detection of C 60 in the PN Ml-1 1 



TABLE 16 

The detected lines in the HDS observations. 



Alab 


Ion 


/(A) /(A) 


Alab 


Ion 


/(A) 


/(A) 


Alab 


Ion 


/(A) /(A) 


(A) 




[7(H/3)=100] 


(A) 






[7(H/3)=100] 


(A) 




[7(H/3)=100] 



3661.21 
3662.26 
3663.40 
3664.68 
3666.10 
3667.71 
3669.46 
3671.38 
3673.74 
3674.84 
3676.36 
3679.35 
3682.81 
3686.83 
3691.55 
3697.15 
3703.65 
3711.97 
3715.08 
3721.94 
3726.03 
3728.81 
3734.37 
3750.15 
3770.63 
3797.90 
3819.60 
3835.38 
3889.05 
3918.97 
3920.68 
3964.73 
3970.07 
4026.18 
4056.91 
4068.60 
4076.35 
4089.30 
4101.73 
4128.66 
4143.76 
4267.15 
4307.44 
4317.19 
4340.46 
4352.52 
4363.21 
4368.24 
4369.56 
4377.77 
4387.93 
4414.90 
4416.98 
4471.47 
4541.59 
4571.10 



H31 

H30 

H29 

H28 

H27 

H26 

H25 

H24 

H23 

Hell 

H22 

H21 

H20 

H19 

H18 

H17 

H16 

H15 

O III 

H14 

[O II] 

[O II] 

H13 

H12 

Hll 

H10 

He I 

H9 

H8 

CII 

C II 

He I 

H7 

He I 

Nil 

[S II] 

[S II] 

O III 

H6 

Nil 

He I 

CII 

Nm? 

C III 

H5 

Nil 

[O III] 

O I 

O in 

Mn II 

He I 

Oil 

Oil 

He I 

Hell 

Me I] 



0.335 
0.335 
0.335 
0.335 
0.334 
0.334 
0.334 
0.333 
0.333 
0.333 
0.332 
0.332 
0.331 
0.330 
0.329 
0.328 
0.327 
0.325 
0.324 
0.323 
0.322 
0.322 
0.321 
0.317 
0.313 
0.307 
0.302 
0.299 
0.286 
0.279 
0.279 
0.267 
0.266 
0.251 
0.243 
0.239 
0.237 
0.233 
0.230 
0.222 
0.217 
0.180 
0.167 
0.164 
0.157 
0.153 
0.149 
0.148 
0.147 
0.145 
0.142 
0.133 
0.132 
0.115 
0.094 
0.084 



0.681 
0.698 
0.483 
0.865 
0.465 
0.583 
1.135 
0.961 
1.182 
0.174 
1.025 
1.384 
1.482 
1.569 
2.015 
1.853 
1.741 
2.328 
0.296 
2.583 

55.380 

19.692 
2.590 
3.438 
3.833 
7.047 
0.363 
8.950 

13.681 
0.242 
0.510 
0.598 

18.792 
0.949 
0.254 
1.790 
0.634 
0.095 

27.822 
0.216 
0.175 
0.614 
0.099 
0.120 

46.900 
0.028 
0.137 
0.317 
0.137 
0.060 
0.238 
0.191 
0.223 
1.718 
0.249 
0.221 



±0.111 
± 0.073 
± 0.062 
± 0.075 
± 0.063 
± 0.079 
± 0.082 
± 0.063 
± 0.075 
± 0.045 
± 0.076 
± 0.074 
± 0.058 
± 0.054 
± 0.097 
± 0.065 
±0.051 
± 0.052 
± 0.041 
± 0.074 
± 0.526 
± 0.302 
± 0.064 
±0.051 
± 0.069 
± 0.065 
± 0.034 
± 0.083 
±0.100 
± 0.032 
± 0.025 
± 0.023 
±0.104 
± 0.020 
± 0.029 
± 0.047 
± 0.055 
± 0.014 
±0.127 
± 0.052 
±0.019 
± 0.027 
±0.019 
±0.018 
±0.174 
±0.011 
±0.015 
± 0.029 
±0.031 
±0.016 
± 0.008 
± 0.017 
± 0.020 
± 0.030 
± 0.021 
± 0.017 



4574.8 

4591.12 

4630.54 

4634.12 

4635.32 

4640.03 

4640.64 

4641.81 

4641.85 

4649.14 

4650.84 

4652.05 

4658.64 

4676.23 

4685.68 

4701.53 

4713.17 

4754.69 

4769.43 

4789.57 

4861.33 

4881.00 

4921.93 

4924.54 

4931.80 

4934.08 

4958.91 

5006.84 

5015.68 

5017.22 

5030.33 

5032.07 

5041.02 

5047.74 

5055.98 

5080.49 

5121.83 

5131.25 

5143.29 

5145.17 

5145.75 

5146.45 

5191.82 

5197.90 

5200.26 

5270.40 

5274.78 

5298.88 

5304.54 

5342.43 

5347.89 

5351.81 

5365.10 

5461.91 

5494.69 

5517.66 



[Mn III] 

Sn 

Nil 

Nm 

Fe II 

C III 

Nm 

Oil 

Nm 

On 

On 

C III 

C IV 

On 
Hen 
[Fe III] 
He I 
[Fe III] 
[Fe III] 
Nil 
H4 

[Fe III] 
He I 
[Fe III] 
[O III] 
Ball 
[O III] 
[O III] 
He I 
Nil 
[Fe IV] 
CII 
Si II 
He I 
Sill 
[Ni III] 
CII 
[Krv] 
[Fe III] 
C II 

[Fe VI]? 
[Fe III] 
[At III] 
[N I] 
[N I] 
[Fe III] 
[Mn II] 
[Fe II] 
C III 
CII 
Oil 
[Cr III] 
[Mn III] 
Fe II 
[Fe III] 
[CI III] 



0.083 
0.078 
0.066 
0.065 
0.065 
0.063 
0.063 
0.063 
0.063 
0.060 
0.060 
0.060 
0.058 
0.053 
0.050 
0.045 
0.042 
0.030 
0.025 
0.017 
0.000 
-0.005 
-0.016 
-0.017 
-0.019 
-0.019 
-0.026 
-0.038 
-0.040 
-0.040 
-0.043 
-0.044 
-0.046 
-0.048 
-0.050 
-0.056 
-0.065 
-0.067 
-0.070 
-0.070 
-0.071 
-0.071 
-0.081 
-0.082 
-0.083 
-0.098 
-0.099 
-0.104 
-0.105 
-0.112 
-0.113 
-0.114 
-0.117 
-0.135 
-0.141 
-0.145 



0.059 
0.082 
0.055 
0.447 
0.033 
0.012 
0.115 
0.026 
0.044 
0.022 
0.031 
0.213 
0.137 
0.096 
0.078 
0.069 
0.260 
0.039 
0.053 
0.051 
100.000 
0.049 
0.464 
0.016 
0.136 
0.016 
6.440 
18.544 
0.989 
0.014 
0.036 
0.106 
0.043 
0.105 
0.119 
0.137 
0.063 
0.140 
0.048 
0.092 
0.018 
0.115 
0.019 
0.201 
0.123 
0.071 
0.055 
0.118 
0.063 
0.032 
0.018 
0.047 
0.022 
0.038 
0.064 
0.029 



± 0.010 
± 0.015 
± 0.007 
± 0.027 
± 0.004 
± 0.001 
± 0.004 
± 0.001 
± 0.002 
± 0.006 
± 0.007 
±0.019 
± 0.008 
± 0.023 
± 0.021 
± 0.008 
± 0.012 
± 0.009 
± 0.006 
± 0.016 
±0.130 
± 0.005 
± 0.007 
± 0.008 
± 0.026 
± 0.004 
±0.031 
± 0.082 
± 0.004 
± 0.004 
± 0.014 
±0.018 
± 0.008 
± 0.006 
±0.011 
± 0.023 
± 0.009 
±0.015 
± 0.003 
± 0.005 
± 0.008 
± 0.007 
± 0.004 
± 0.009 
± 0.008 
± 0.004 
± 0.009 
± 0.009 
±0.016 
± 0.004 
± 0.005 
± 0.006 
± 0.004 
± 0.021 
± 0.014 
± 0.006 



5537.60 
5554.83 
5577.20 
5659.60 
5666.63 
5690.43 
5695.92 
5730.66 
5754.60 
5826.42 
5831.70 
5875.66 
5958.39 
6046.23 
6074.20 
6077.90 
6267.81 
6300.34 
6312.10 
6347.09 
6363.78 
6527.24 
6548.10 
6562.77 
6578.01 
6583.50 
6678.16 
6688.79 
6701.49 
6716.44 
6730.82 
6744.33 
7002.17 
7018.63 
7037.25 
7065.18 
7065.71 
7101.07 
7135.80 
7179.88 
7231.33 
7231.48 
7236.20 
7236.42 
7236.80 
7237.17 
7254.53 
7281.35 
7318.92 
7319.99 
7329.67 
7330.73 
7442.20 
7452.60 
7457.90 
7468.31 



[CI III] 

Ol 

[O I] 

Cm 

Nil 

Si I 

C III 

Nil 

[Nil] 

C III 

Nm? 

He I 

Ol 

Ol 

CI 

CII 

[V II]? 

[O I] 

[S III] 

Si II 

[O I] 

[Nil] 

[Nil] 

H3 

CII 

[Nil] 

He I 

Cl 

Nil 

[S II] 

[S II] 

Cm 

Oi 

Ca I]? 

C III 

He I 

He I 

Cal 

[Ar III] 

Ca I 

CII 

Cal 

Cal 

CII 

CI] 

CII 

Ol 

He I 

[Oil] 

[On] 

[On] 

[On] 

C IV 
[Fe II] 
CII] 
Ni 



3.149 
-0.152 
-0.156 
-0.169 
-0.171 
-0.174 
-0.175 
-0.181 
-0.185 
-0.196 
-0.196 
-0.203 
-0.215 
-0.228 
-0.232 
-0.232 
-0.258 
-0.263 
-0.264 
-0.269 
-0.271 
-0.293 
-0.296 
-0.298 
-0.300 
-0.300 
-0.313 
-0.314 
-0.316 
-0.318 
-0.320 
-0.322 
-0.356 
-0.358 
-0.361 
-0.364 
-0.364 
-0.369 
-0.374 
-0.380 
-0.387 
-0.387 
-0.387 
-0.387 
-0.387 
-0.387 
-0.390 
-0.393 
-0.398 
-0.398 
-0.400 
-0.400 
-0.415 
-0.416 
-0.417 
-0.418 



0.072 
0.100 
0.040 
0.054 
0.101 
0.038 
2.113 
0.043 
5.986 
0.122 
0.041 
6.155 
0.129 
0.239 
0.057 
0.075 
0.041 
2.203 
0.323 
0.126 
0.748 
0.030 
58.851 

285.000 
0.499 

190.199 
1.600 
0.050 
0.067 
0.442 
1.084 
0.125 
0.210 
0.023 
0.139 
2.276 
0.334 
0.020 
2.457 
0.066 
0.216 
0.077 
0.058 
0.436 
0.057 
0.071 
0.263 
0.319 
3.837 
11.165 
5.513 
5.284 
0.076 
0.119 
0.044 
0.103 



± 0.012 
± 0.023 
± 0.008 
±0.011 
± 0.048 
± 0.017 
± 0.036 
± 0.013 
± 0.077 
± 0.024 
± 0.012 
± 0.083 
± 0.021 
± 0.028 
± 0.014 
± 0.021 
± 0.010 
± 0.041 
± 0.015 
± 0.017 
± 0.021 
± 0.004 
± 0.989 
±4.160 
± 0.012 
± 4.767 
± 0.028 
± 0.012 
± 0.008 
± 0.016 
± 0.027 
± 0.017 
± 0.009 
± 0.006 
± 0.013 
± 0.036 
±0.011 
± 0.006 
± 0.042 
± 0.006 
± 0.012 
± 0.007 
± 0.004 
±0.011 
± 0.004 
± 0.005 
± 0.012 
± 0.013 
± 0.077 
±0.193 
±0.103 
±0.102 
± 0.007 
± 0.010 
± 0.005 
± 0.007 
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TABLE 17 

The detected lines in the ISLE observations. 



A vac 


Ion 


/(A) 


F(X) 


/(A) 


(fxm) 






(erg s 1 crrT 2 ) 


(/(Pa/3,Br7)=100) 


1.129 


O I 


-0.715 


"7 1 A / 1 A \ I A C 1 / 1 CA 

7.34(-14) ± 4.51(-15) 


OCTC I A 10 

2.535 ± 0.188 


1.147 


[P II] 


-0.722 


2.23(-14) ± 3.20(-15) 


a Trr I A110 

0.755 ±0.113 


1 1 AO 

1. 102 


IT o A C { 1 \ 

H2 2-U o(l ) 


—0.728 


ci / 1 c\ _i_ n oa/ 1 A\ 
o.Dl(— ID) ± y.oU(— 10) 


A OOQ _L /I AQC 


1.175 


Ol? 


-0.733 


1 If! 1 /I \ 1 1 O C\ / 1 ^\ 

l.76(-14) ± l.29(-15) 


0.577 ± 0.049 


1 1 OA 

1.1 00 


it 1 1 C/IA 

H2 3-1 o(3) 


-0.736 


1 at/ 1 a \ _i_ n aa/ 1 A\ 
1.02(-14) ± y.OO(-lO) 


A COC _L /I /1QO 

U.D2D ± U.U3o 


1.189 


[Pll] 


-0.737 


s cn/ 1 /I \ _l_ 1 AA/ 1 C\ 

0.50(-14) ± 1.00(-15) 


O 1 A/C I A ("1AC 

2.106 ± 0.095 


1.190 


TT 1 (1 O /A\ 

H2 2-0 S(U) 


-0.738 


A /I/;/ 1 C\ _|_ 1 AA/ 1 c\ 

0.40(-15) ± 1.00(-15) 


a inn 1 ci at c 

0.20y ± 0.035 


1 1 (\A 

i . iy4 


L.al : 


n Tin 
-0. /39 


r n/ 1 C\ _L 1 111 1 CA 

D.12(— ID) ± 1.3 /(-ID) 


A 1 AC J_ <\ <\AZ 

U.10D ± U.U4D 


i.iy / 


U,, 1 

rie I 


A 7A{\ 
—0. /4U 


no/ 1 c\ _i_ A(\i 1 a\ 
o.UV(— ID) ± o.4U(— 10) 


a o/;n J_ A ATA 
U.20U ± U.U2y 


1.199 


N III? 


-0.741 


7.32(-15) ± O.oO(-lO) 


A 1 Q C _L A AT /I 

0.235 ± 0.024 


1 one 
1.2U5 


H 2 3-1 S(2) 


A HA A 

-0. /44 


q Ait i c\ 1 7 nnc 1 
S.43(-1D) ± /.yu(-io) 


A OAO _L /^ /iTO 


1.219 


C I? 


-0.748 


4.09(-lD) ± 7.90(-10) 


a 1 An _i_ a m<; 
0.147 ± 0.020 


1.221 


Ca I? 


-0.748 


A A f 1 C\ I 1 C\f\t 1 /; \ 

3.44(-lD) ± 7.y0(-10) 


A 1 AO I A AO C 

0.108 ± 0.025 


1.226 


TT A 1 C 1 / C \ 

H2 4-2 S(5) 


-0.750 


4.79(-15) ± 6.80(-16) 


a 1 ca i f\ 

0.150 ± 0.022 


Y.ZZy 


L. II 


-0. /51 


1. UD(— 14) ± 0.0U(— 10) 


A Q17 _L /I AOC 
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Si I? 


-0.751 


7.11(— 15) ± 0.00(-10) 


A ti 1 I a nn 
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1.233 


TT 1 IC/1\ 

H2 3-lo(l) 


-0.752 


2.07(-14) ± o.o0(-10) 


A 1 A Al A 

0.042 ± 0.03y 
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T T OA / 1 \ 

H 2 2-0 Q(l) 


-0.754 
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1.5U(— 14) ± 1.08(-15) 
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1 .242 
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H2 4-2 S(4) 
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/.44(-lD) ± l.UD(— ID) 


A TOO J_ A AQ/l 




■ TT '"I n/~\/ o \ 

+ H2 2-0(J(2) 








1.247 


Si I? 


-0.757 
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4.84(-15) ± 2.08(-15) 
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1.24/ 
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He I 


-0.891 
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He I 
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2.122 


H2 1-0 S(l) 


-0.897 


^ nn/ 1 _i_ 1 nc/ 1 c\ 
4.90(-14) ± 1.05(-15) 


O /I "7 1 _L A A AO 

8.471 ± 0.408 


2.138 


Mg II 


-0.898 


4.67(-15) ± 5.30(-16) 


A OAC 1 A 1 A A 

0.805 ±0.100 


2.144 


Mgll 


-0.898 


4.88(-15) ± 8.40(-16) 
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2.154 


H 2 2-1 S(2) 


-0.899 


8.00(-15)±8.10(-16) 


1.373 ±0.156 


2.162 


He 1 


-0.900 


9.78(-15) ± 3.30(-16) 


1.676 ±0.103 


2.166 


H I (Br7) 


-0.900 


5.84(-13) ± 5.23(-15) 
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2.171 


Ol 


-0.900 


8.54(-15) ± 1.08(-15) 


1.461 ±0.199 


2.189 


He II 


-0.902 


6.50(-15) ± 7.60(-16) 


1.108 ±0.141 


2.199 


[Kr III] 
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3.215 ± 0.265 


2.201 


H 2 3-2 S(3) 
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0.669 ± 0.303 
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